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I. INTRODUCTION 



On a June night in 1940, a lone ANSON aircna-ft -flew a 
solitary path through the skies o-f East Anglia, England. The 
wireless operator was listening -for traces o-f signals 
propagating through the dark space around the aircra-ft. Near 
an important military target, a signal came through loud and 
clear in his earphones. First, a series of dots at the rate 
of sixty per minute. Then, as the ANStS^I droned onward, the 
dots disappeared into a steady tone. A little later, the 
steady tone broke up into a series of dashes, at the rate of 
sixty per minute. 

This signal , known as a Lorentz beam, was transmitted from 
Germany and it posed a considerable threat to the high value 
military target. The signal was really composed of two beams, 
one beam of dots, and one beam of dashes, shown in Figure 1.1. 
By using the reception of the dots and dashes as an indication 
of whether they were on course, German bombers could navigate 
to their targets. The phasing of the dots and dashes was such 
that they combined into one steady tone in the overlap. 

Receipt of the steady tone meant the aircraft was on the 
course laid dcx«jn by the transmitting site and receipt of dots 
or dashes meant the bomber was right or left of the beam 
overlap. With that information, the bomber crews could easily 
get back on course and perform their mission. C13 
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The beams were ultimately countered by the British using a 
variety o-f methods, but the invaluable intelligence gained by 
the ANSON bomber over East Anglia was only the beginning o-f a 
wartime e-f-fort that was to save thousands o-f lives, both 
civilian and military. That brave and skillful act has been 
repeated hundreds of times since, in a "police action" and an 
undeclared war, but the priceless results have been the same, 
that of saving lives. 



Kzf- Ste3(fy note zone 




Dotj and dashes interlock 
so that Mfhere beams overlap 

a steady note is heard THE LOHENZ BEAM 

Figure 1.1. The Lorentz Beam 

Since the arrival of the first electronic transmitter and 
receiver, man has been listening to his own radiations for 
communications, pleasure, and information. World War II, 
however, brought a new maturity to listening to his own 
radiations. Now he was listening to an opponent's 
transmissions to derive and use information about them. 
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About 1969, the term Electronic Support Measures < ESM) was 
coined C2; p.571 and now reflects the important business of 
obtaining and using in-formation gained on the electromagnetic 
radiations o-f another. The alert operator aboard the lone 
British aircra-ft used his early and unde-fined ESM receiver to 
obtain knowledge that was used against the German 
transmi t ters . 

A. THE ELECTRONIC WARFARE <EW) ENUIRONMENT 

0-f all the sensors that an enemy might use, the most 
important, at least as an electronic threat, is radar. This 
sensor is almost always the long and short range eyes and ears 
o-f the opposing -force, and denial of its use gives a strategic 
and tactical advantage to the person capable of doing so. C2: 
p.7] 

Consider the case of a penetrating bomber. In order to 
jam, deceive, or otherwise negate the effectiveness of the 
radars of the defense system, emissions must be detected, 
analyzed, and classified as to their threat. In this 
endeavor, the ESM receivers being used will not lack for 
electromagnetic sources to perform their analyses. Early 
warning, target detection, acquisition, and tracking, or even 
missile and interceptor guidance and control radars will paint 
the aircraft with an invisible illumination. 

These electromagnetic energies will be of many 
frequencies, antenna scans, pulse shapes, and modulations, and 
even though radar systems are well understood and sometimes 
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easily countered, the sheer numbers and types ohF received 
signals will tax the capabilities o-f the ESM system. 

As the bomber nears the target, the signal densities will 
soar because the de-fense is toughest close to the target. For 
instance, it has been reported that over 108 SAM sites were 
once located on the Egyptian side of the Suez Canal. As each 
battery transmits several signals, easily 390-509 signals 
could be found from this one class of weapon. When 
communication and radio location signals are added, the signal 
density problem is even further complicated. t2: pp. 27-333 

Quite easily, the combat EW environment that the bomber is 
in will reach 500,900 pulses per second. [3; p.543 

Figure 1.2 C43 gives an indication of the pulse densities 
the bomber might encounter. The pulse density is dependent 
upon both system sensitivity and platform altitude but is not 
uniformly dependent on these parameters. At low sensitivities 
<-20 to -70 dBm) pulse density varies mostly with sensitivity 
and not with altitude. At high sensitivity 0-100 dBm) pulse 
density is almost totally altitude dependent. 

The identification of various emitters and the evaluation 
of the possible threats they pose to the aircraft requires the 
sorting of these pulses. This separating or deinterleaving 
of signals is a difficult problem that is discussed in Chapter 
II . 
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PULSE DENSITY 




Figure 1.2. Pulse Density vs. Sensitivity and Altitude 
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B. ELECTRONIC WARFARE SUPPORT MEASURES <ESM) 



The broad term Electronic War-fare (EW), or sometimes 
re-ferred to as Electromagnetic War-fare <EMW) , is used to 
denote military actions that involve the -friendly exploitation 
o-f the electromagnetic spectrum and the prevention o-f hostile 
use o-f the el ec tromaget i c spectrum. It should be realized 
that the electromagnetic spectrum covers -from zero -frequency 
to the in-finite -frequency, which includes the optical, 
in-frared, and laser frequencies. C53 Coverage of these 
additions to the RF spectrum adds considerably to the problem. 

EW may be divided into three major divisions: 

1. Electronic Warfare Support Measures <ESM) 

2. Electronic Counter Measures <ECM) 

3. Electronic Counter Counter Measures (ECCM) 

Electronic Support Measures, the usual shortened form of 
Electronic Warfare Support Measures, will be defined shortly. 
ECM, quite briefly, is the generation of intentional 
electronic interference between electronic systems for a 
useful purpose. ECCM , also quite briefly, is action taken to 
counter the detrimental effects of the opponent''s ECM. I6t 
p. 13 

This section will deal only with the ESM arena. 

1 . A Definition of E^ 

Consider a specific definition of E^i. ESM is; 

That division of Electronics Warfare involving actions to 
search for, intercept, locate, record, and analyze radiated 
electromagnetic energy for the purpose of exploiting such 
radiations in support of military operations. C7; p.63 
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This definition carries an underlying statement: ESM 
is the source of information required to carry out the other 
Electronic Warfare divisions as well as threat detection, 
warning, avoidance, target acquisition, and homing. C4: 
pp.3“<63 In addition, by examining the individual actions 
specified in the definition we can gain further insight into 
the ESM arena. 

The search for and interception of electromagnetic 
radiations is simply a basic part of the ESM mission. The 
intercepted signal must be analyzed and a determination made 
of its character, parameters, and location. Locating the 
emitter is usually a function for direction finding antennas. 
Recording of the signal, certainly the parameters, is usually 
done simultaneously with the analysis. On some missions a 
recording is made of the signal for postflight analysis with a 
mimimum of information being displayed to the operator in real 
time. Generally it is desired to produce this information 
from an ESM system with an inherent automatic processing 
capability that detects, classifies, and flashes a warning if 
need be. C2: p.36] 

The information that the ESM mission provides ranges 
from that collected with regularly scheduled and dedicated 
collection excursions to debriefs of crews after completion of 
a routine patrol mission. To a tactical commander, this 
information is more extensive and useful than plain 
intelligence. It provides the precise state of the electronic 
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defense, including technical characteristics and emitter 
location, and becomes part of mission preflight briefings. 

In obvious fact, radars not transmitting are radars 
that can't be intercepted and analyzed. Thus, if there are no 
transmissions from specific desired radars, or any at all, 
tactics that might provoke them into radiating are viable in 
the ESM scenario. Indeed, this action is included in the 
basic rules for Electronic Reconnaissance. C2: p.651 
2 . The Role of E^ in Electronic Warfare 

In general, ESM is passive electronic warfare, ECM is 
active electronic warfare, and ECCM may be either. These 
distinctions may not be clear cut, however. ESM may involve 
the radiation of a signal to determine the char ac ter i st i cs or 
electronic reaction of enemy equipment. ECM may involve the 
passive reception of enemy signals in order to decide which 
one to counter C2; p.33. The latter could be the use of a 
tactical E^ system integrated with an ECM system for platform 
protection. The relationships between E^l, ECM, and ECCM are 
shown in Figure 1.3. C2: p.23 

Enemy Combat ’ Friendly Reconnaissance 




Figure 1.3. The Interactions of EW 
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The previous consideration o-f the penetrating bomber 
can be continued in the context o-f ESM"s role in electronic 
war-fare . 

The bomber will almost surely face an integrated 
defense network. It is absolutely essential to know 
beforehand the technical char ac ter i st i cs and location of the 
electronic defense that will be used against the bomber. Some 
kind of electronic reconnaissance information is needed in 
order to maximize the effectiveness of ECM and ECCM equipment. 
Electronic intelligence (ELINT) and ESM are our sources of 
informat i on . 

The value of the knowledge gained from electronic 
reconnaissance of the defense system is proportional to its 
currency. The use of information from electronic intelligence 
which was gained from many sources and over a longer period of 
time is important, but the immediate needs of the tactical 
commander are met using ESM. This information is collected in 
the hours before takeoff, undergoes minimal analysis, and is 
used for mission planning. C2: p.373 Thus, ESM is electronic 
reconnaissance to determine the present state of the defense 
system . 

ELINT, which has been subjected to extensive analysis, 
has been used to design and develop the ECM suite of the 
bomber. But the tactical commander still uses ESM to plan and 
accomplish the mission. Though they may share a corronon 
collection platform, the difference between them can be seen 
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in their operational control , end use, and time or duration of 
collection. ESM is used by the tactical commander today, and 
ELINT goes to the rear echelon for further study. l 2: p.74] 
Figure 1.4 C8] shows the operational role of ESM in 
electronic warfare. Note the central position of ESM as it 
"feeds" both ECM and ECCM and operational readiness. 
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C . SUMMARY 



The previous sections only touch upon the vast nature of 
EW and ESM. They provide a general groundwork -for 
understanding the purpose of using ESM and its role in the 
electromagnetic conflict. 

This conflict is extremely important. Nei/* and modern 
weapon systems depend upon victory in the electromagnetic 
conflict as a prerequisite for victory in battle. C2; p.lD 

Underscoring this premise are new policies and procedures 
by the Defense Department which will accelerate development 
and introduction of electronic warfare systems. The rapid 
escalation and introduction of new Soviet radars places 
current EW sytems at a disadvantage, especially in countering 
fire control and missile guidance radars. C9] 

The future probably holds a greater integration of ESM 
(certainly tactical ESM) and ECM systems, with design for 
military aircraft leading the way. Already systems that 
combine EDi capability with passive warning equipment for 
aircraft are being sought by many NATO nations, such as the 
British Royal Air Force C103, the Belgian Air Force C113, as 
well as the United States Air Force [121. In addition, and 
with considerable impact on system design, a trend toward 
carrying these systems internally rather than in pods hung 
underneath the aircraft seems to be evident C 18 and 121. 

In the hardware area, the use of embedded computers and 
microprocessors with high use of ^v*HSIC technology will become 
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prevalent, along with all o-f the rami -f i cat i ons o-f computer 
use. The sophisticated signal processing needs o-f B'l will have 
to be met, with veri-fiable, reliable, maintainable so-ftware, 
and the development o-f high capacity, non-volatile data 
storage . C 13] 

Upcoming and -foreseeable developments are EW systems -for 
space applications C2: pp. 171-178]. All o-f the unique 
environmental conditions o-f space will generate new and 
di-f-ferent challenges -for ECM and ECCM. These must be met in 
order -for manned space stations and satellites to survive. 

D. SCOPE OF THESIS 

A-fter these introductory sections Chapter II examines some 
o-f the system design requirements of an ESM system. The 
various components are looked at in some detail. Signal 
processing is addressed, beginning with basic parameter 
measurement, and ending with a discussion of our basic 
problem, that of deinterleaving received pulse trains. 

The attack on the problem begins in Chapter III, with a 
discussion of the theory of Walsh and Rademacher Functions. 
Application of this theory is carried out in Chapter IV, where 
the Walsh transform is used on the deinterleaving problem. 
Finally, conclusions are developed, and listings of the 
programs used are given in order that follow on work can be 
done . 

The attack on the problem using the Walsh Transform will 
begin with basic examples and proceed with an application of 
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the trans-form directly to an interleaved pulse train. The 
goal is to determine whether the Walsh Trans-form outputs a 
unique -feature that can be used in a deinterleaving algorithm. 

The radar pulses can be represented by time o-f arrival 
(TOA) tags, which are output by the ESM receiver to the 
preprocessor. A TOA tag can be thought o-f as merely a “mark" 
in time that indicates the receipt of a pulse. The stream of 
incoming pulses forms a string of TOA TAG'S. The Transform is 
applied to the TOA tag string. 

The Walsh Transform is easily applied to binary or two 
level functions,- and some promise is held for its use in a 
deinterleaving algorithm if a PRI/PRF recognition feature 
exists in the transform. 

With the speed and ease of computation of the Walsh 
Transform, a processor based on these functions should be fast 
and simple. Consider the programs used to compute the Walsh 
transforms in this thesis. They were originally written in 
FORTRAN, and were easily converted to BASIC for use in this 
thesis. 

This simple, easy to use language normally brings a 
response that would question the use of such an “elementary" 
language. That response is misinformed. 

BASIC is easy to learn, and easy to use. It performs well 
when compared to other structured languages, and has excellent 
mathematical and scientific processing capability. The use of 
the language certainly has not hindered the growth of the 
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computing power available in a mi crocompu ter . The author''s 
personal computer, the IBM Personal Computer, with 64K of RAM, 
has the same amount of memory that a -frontline U.S. Navy P3-C 
Orion ASW aircra-ft uses to process seven tactical work 
stations, plus handling all the navigation, displays, and 
other systems. 

Throughout this thesis, a microcomputer was used to 
compute all the results given, and only the lack o-f a thesis 
quality graphing device prevented the author -from doing all 
work on his own personal computer. Don't underestimate a 
microcomputer's capabilities. The day is coming when 
main-frame-like computing power will be available in each 
o-f-fice as close as the desk or table. 
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II. ELECTRONIC SUPPORT MEASURES SYSTEMS 



The ESM system used on the modern battlefield must be able 
to give rapid and accurate assessments of the complex RF 
environment. Detection and recognition of hostile radar 
transmissions must take only a few seconds because to do less 
invites destruction or the loss of the ability to deal a 
decisive blow to the enemy. C3; p.54I 

The design of ESM systems, in view of their importance, 
would seem to follow a rational and orderly path from concept 
to development to operational capability. 

This usually isn't the case, unfortunately, as the 
reactive nature of EW will impede an orderly flow or 
development. Without perfect knowledge of an enemy's radar 
and EW capabilities, surprises occur which tend to drive the 
normal development into a "reaction" driven process. The 
normal process also tends to be delayed by a desire to wait 
until a precise definition of the threat is in hand. By the 
time it is, a reaction process is needed to counter the 
threat. £2; p.533 

As referred to in the Introduction, page 25, this design 
process is undergoing changes to improve the response and 
lessen the reaction design. Changes in the design process 
will not change the design requirements. 
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The system is designed to intercept many di-f-ferent 

electromagnetic signals. Ideally, the system should be able 
to accomplish as quickly as possible the -following tasks: 

1. Intercept a transmitted signal at any -frequency 

2. Determine the type o-f modulation in the signal 

3. Identi-fy the usable intelligence carried in the signal 

4. Measure the direction o-f arrival o-f the wave-form -for 
calculations that locate the transmitter. 

5. Process and record the signal characteristics 

6 . Display to the operator and trans-fer to the computer -for 
decision making 

In short, the system must gather, process, and display all 
signals o-f interest. 12: p.l913 

Meeting these requirements requires a lot o-f attention and 
hardware. No single tuner can cover the area o-f interest 
(usually 500 MHz to 30 GHz), so many tuners, antennas, and 
other system units are required. 

With this equipment, an ESM system approaches the problem 
with these basic intents: 

1. Intercept the signal . Conversion to a usable -form must 
now take place immediately. An alarm should be given to 
indicate a signal has been intercepted. Some analysis 
or sorting must be done immediately to determine i-f the 
signal is o-f high threat, since signal -frequency and 
modulation usually correlate with degree o-f threat. 12; 
p. 1913 

2. Analyze the signal -for -further in-formation. Once it's 
decided there isn't an immediate threat, other 
parameters o-f the signal can be determined -for recording 
and later decision making. 

Each o-f these approaches has more or less importance 
depending upon the plat-forms. An airborne plat-form has less 
space to accommodate equipment -for reception and analysis, but 
is usually exposed to more immediate threats. There-fore, an 
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airborne plat-form devotes more attention to interception and 
iiwnediate analysis than the ground plat-form, which probably 
is more concerned with a -full analysis and recording. 
Hope-fully, this entire thesis has the flavor of applying to an 
airborne platform. 

An examination of a typical airborne ESM system will begin 
in the next section. 

A . COMPONENT SYSTEMS 

Figure 2.1 C43 shows a typical ESM system. Of particular 
interest and discussion in this chapter will be the receiving 
system, the preprocessor , and the processor. The Data Files 
(Active Emitter and Emitter Parameters), are used by the 
system as records for current signals being received (Active 
Emitter) and comparison of parameters from processed signals 
with previously entered (a priori) parameters of signals of 
interest or expected signals (Emitter Parameters) . Electronic 
Counter-Measures (ECM) handoff and equipxnent is not really a 
part of the E3i system, but is included as a function that 
could be available for manual/automatic Jamming or deception 
of high threat signals. 

1 . The Receiving System 

Inseparable in their functions, the antennas and the 
receiver form the interception and detection system that is at 
the forefront of the entire process. This section will cover 
a few of the important aspects of the antennas and receiver 
systems. 
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Figure 2.1. An ESM System. 
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a. Antennas 



The energy received by any plat-form is a very 
small -fraction o-f that transmitted. Consider Table 2.1 C14: 
p. 9-331. Notice the numbers in the received power column are 
magnitudes below one milliwatt in almost all instances. 
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High sensitivity in the system is needed to use 
this power. Good antenna design is the starting point -for 
reception o-f the signals that are desired to be analyzed, and 
good design is very dependent upon the application. 

The basic equation that is -fundamental to the 
amount o-f power available -for use is; 

Pp = Pt6^.x2/<<410 2R2) <2.1) 

where the variables and units are as indicated below: 
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Pp = received signal power in watts 
= transmitted power in watts 
= e-f-fective transmitting antenna gain 
Gp = effective receiving antenna gain 
X = wavelength (same units as R) 

R = separation between transmitter and receiver 

In this age of "if you can see it, you can hit it" 
accuracy with guided missiles, it is not recommended that 
shorter ranges be used to improve the amount of power being 
received. Thus, the only other available factor under the 
designer's control (besides receiver design) is the receiving 
antenna gain. 

Our antenna gain is affected by several factors. 

G = 4'SAP/X2 (2.2) 

where A and P are defined 

A = capture area 
P = antenna efficiency 

The capture area is the physical size of the 
antenna. The efficiency is an inherent property and a 
function of the type of antenna (usually between 0.5 and 
0 . 6 ) . 

The beamwidth of the pattern that has the 
designer's interest also has an effect on the antenna gain 
C 15] . It is an inverse relationship C 14; p.29-3] that 
depends upon the current distribution across the aperture. 

A typical reflector antenna, for instance, has 
G ss 20 ,000/(eg#0) (2.3) 

where 00 and ^0 are the ha 1 f -power beamwidths in degrees 
measured in the principal planes of the pattern. 



34 



The Gain 6 is the power gain, and should be used 
in radar equations because it includes losses introduced by 
the antenna (Eq. 2.2). Bearing resolution and coverage 
•functions -for direction -finding would lead to a consideration 
o-f the directive gain. The directive gain is more descriptive 
o-f the antenna pattern. 

A use-ful expression -for directional gain is 

Gq = 4K/<ee*B^ <2. 4) 

where 63 and ♦g are the hal -f -power beamwi ths in radians. 

[ 153 

Gain must be considered in terms o-f signal 
reception. Two o-f the system's tasks are to receive the 
signals which may come -from any direction (ornni 
considerations) and locate the emitter <DF considerations). 

The two tasks are not necessarily complementary, so a 
trade-o-f-f in system per-formance can be accepted or 2 antennas 
can be used, one -for each mission. 

For an airborne plat-form, weight, space, and 
structural requirements may limit the size o-f an ESM antenna, 
but the altitude o-f the aircra-ft will usually o-f-fset some o-f 
these disadvantages. Dipoles, slotted monopoles, and sur-face 
wave antennas C14; p. 29— 303 are antennas that can be used in 
airborne applications. 

The DF -function o-f the system could use any o-f the 
mentioned types o-f antennas i-f the system uses amplitude or 
phase comparisons between two collectors. [14: p. 29-3 13 
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Loop antennas are a common amplitude sensitive DF 
antenna o-ften used in a system called a goniwneter [16: p. 
2h-14]. Horizontal dipoles or dipoles combined with 
re-flectors are also antenna types that are used in DF systems, 
as well as Adcock arrays C 16: p.2h-18]. 

Generally, it is desired to have high gain in 
antennas -for an system, and the speci-fic application will 

dictate the type and size o-f antenna that can be used. More 
o-ften than not, an omni-directional antenna o-f some type is 
used -for search, and a rotating, -fixed, or multi-port highly 
directional antenna- is used -for direction -finding, 
b. Receivers 

Most discussions on EW wi 1 1 address the receiver 
technology used to detect the signal brought in by the antenna 
system. Several types of receivers are in use by the Air 
Force, Navy, and Army, with the mission usually dictating the 
choice of technology. 

A quick and general overview of receiver types and 
technologies follows in this section. Although some of the 
technologies have been around for some time, and some still 
yet to come, they are increasingly being dominated by the 
processing requirements and the impact of microcomputers. 

Cost saving and simplification of the microwave design 
requirements are a result of this increasing influence. In 
addition, combinations of receiver technologies are being used 
to overcome limitations of individual receiver concepts. 
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The future holds new technologies for the detection of 
signals and increases of speed and bandwidth of single 
receiving devices. Surface Acoustic Wave devices will 
certainly play an important role, as well as acousto-optic 
techniques. Increases in computational speed and capacity may 
aid or deter some technologies by eliminating some of the need 
for higher sensitivity in the detection scheme. 

Table 2.2 C 17] gives some typical ESM receiver 
spec i f i cat i ons . 



Table 2.2. Typical Receiver Specifications. 



Frequency Range 
Signal type 
Sensi t ivi ty 
Resol u t i on 

Amplitude accuracy 
Bearing accuracy 
Pulse width resolution 
Time of Arrival resolution 
Signal density 
Intercept probability 



0.5 to 18 GHz 
100 ns pulse to CW 
<-70 dBm 
5 MHz < pu 1 se) 

1 MHz (CW) 

1 dB 
5 - RMS 
100 ns 
50 ns 
106 pps 
1007. 



(1) Direct Detection . Simplicity in design gives 
the direct detection receiver advantages in cost, reliability, 
size, and weight. This design is also proven, currently in use 
in the front end of radar warning receivers <ra*lR) , giving a 
high Probability of Intercept (POD for signals above its 
Minimum Discernable Signal. 

Today's high density environments pose 
problems for the DD receivers, degrading their performance 
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quicKly. Frequency resolution is not very high as the 
bandwidth is usually an octave or more in width. [18; p.262 
Sensitivity is a problem, although RF amplification can 
provide improvements. 

Quite susceptible to ECCM, the DD receiver 
cannot handle PRI agile signals, and does not provide 
frequency measurement, an important sorting parameter. 

Figure 2.2 is a direct detection receiver. 




Figure 2.2. Direct Detection Receiver 



< 2) Instantaneous Frequency Measurement (IFli) . A 
form of a direct detection crystal video receiver, the IFM 
provides frequency discrimination and measurement with a high 
POI . Often an output of the receiver is a polar display of 
amplitude versus frequency C16: p.2a-123 and the term Polar 
Discriminator is used in place of IFM. 
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The operation o-f the discriminator is the 
heart o-f the receiver. Note Figure 2.3(a). Input signals 
are divided into two parts and phase correlated a-fter one 
signal is delayed by a known time. The length o-f the delay 
line is related to the amount of frequency resolution desired 
and the wavelength that corresponds to the desired bandwidth. 
By banking several di scr imi nators and different lengths of 
delay line, as shown in Figure 2.3(b), a desired frequency 
band can be covered. This approach allows the use of 
practical and relatively inexpensive delays rather than 
attempting to use one delay line precise enough to handle the 
desired frequency region. The length of the longest delay 
becomes the frequency resolution of the bank. CIS: pp. 32-333 

Measurement of the frequency of the inccxning 
signal with these analog discriminators (or current technology 
Digital IFNI'S) allows this important parameter to simplify the 
sorting process. 

On the down side, the IFM doesn't handle high 
density environments very well unless a subsystem that 
determines when 2 pulses occur simultaneously is included. 
Also, the delay lines require special attention (constant 
temperature) and this usually increases the required input 
power to the system. 

Working in conjunction with an analysis 
receiver , the IFM makes an ideal acquisition receiver by 
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(a) 




(b) 



Figure 2.3. IFM Receiver Concepts. 
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providing the -frequency of a signal of interest. Analysis 
receivers of another type can then be tuned to the frequency. 
C18; p. 32-343 

(3) Superheterodyne . Most common of the 
receivers in current use, the superheterodyne provides high 
selectivity, good resistance to jamming, and is a proven 
desi gn . 

Drawbacks include slow time in scanning the 
frequency band of interest and inability to see frequency 
agile signals. C19; p.4593 




Figure 2.6. Superheterodyne Receiver. 

<4) Channel i zed . Efforts to increase the PCI of 
superheterodyne receivers and deal with an inability to handle 
high density environments led to the concept of a bank of 
filters followed by individual signal detectors that 
determine when a signal is within the filter bandwidth. The 
RF spectrum is simply broken into pieces. See Figure 2.5. 
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A key advantage of the channelized receiver is 
removal of pulse overlap of near simultaneous signals <only 
pulse frequencies that exceed the channel spacing) . C17: 

p. 185] 

By following an initial bank of filters with 
multiplexers and fixed oscillators, then repeating the 
scheme with a smaller bandwidth filter, the RF input can be 
down converted to a final baseband frequency range with a 
higher frequency resolution. 




Figure 2.5. Channelized Receiver 

<5) Ccffnpressive . This receiver can be 
character ized as a fast scanning superheterodyne receiver. 
Ideally, the local oscillator scans the RF bandwidth being 
covered in a time less than the narrowest pulse to be 
intercepted. Its high POI and the ability to handle wideband 
signals and frequency agile signals makes this receiver a good 
choice for an ESM receiving system. It has an excellent 
ability to separate signals closely spaced in frequency. 
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B&sically, the compressive receiver provides 
an IF signal that is up-chirped in -frequency -from the RF 
input. The compressive -filter operates in such a way that the 
low -frequency components o-f the IF signal are delayed longer 
than the high -frequency components, with the result that the 
shape o-f the output pulse of the detector is in essence the 
Fourier Transform of the IF signal pulse. The maximum 
amplitude of the detector output falls at a point in the scan 
time that is proportional to the frequency of the RF input. 
CIS; pp. 30-31] 

Criticality of the alignment between sweep and 
compressive delay makes this complex technique difficult to 
manage C 19; p. 459]. But it is at the forefront of some of 
the newest technologies, with SAW's performing the delay 
necessary for the up and down chirp. Note Figure 2.6. 
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Figure 2.6. Compressive Receiver. 



<6) Bragg Cell . An acousto-optic Bragg cell will 
interact a propagating acoustic wave with an optical beam to 
give a diffracted output proportional to the frequencies and 
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power present in the RF input signal. The acoustic wave, 
which varies the index o-f re-fraction o-f the cell material, is 
generated by applying the RF input to an acoustic transducer 
in the Bragg cell. A laser beam is de-flected and modulated in 
intensity by the index variations which are proportional to 
the RF input -frequency. A Fourier Trans-form lens collects 
the proportionately di-f-fracted light which then -falls on a 
photodetector array. The output o-f the detector array is 
similarly proportional to the RF input -frequency. C 18; p.283 
Having pluses -for high sensitivity, high POI , and 
high selectivity makes this receiver a choice o-f the -future. 

It doesn't handle -frequency agile sources well, however, 
because o-f the time needed in determining the output o-f the 
photodetector array. Pulse width measurement is also 
eliminated with the use o-f the detector array because the 
array only responds to the energy contained in several pulses 
-from the lens. C19; p.A621 Figure 2.7 is the receiver . 




Figure 2.7. The Bragg Cell Receiver. 
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2. The Processing System 



A steady stream o-f pulse descriptor words from the 
receiving system must be analyzed for the vital information 
that is contained in them. The processing system must extract 
this information and make decisions regarding the emitters 
whose signals it has processed. Information extraction, 
analysis and identification must be fast to avoid a mission 
and life threatening situation. 

This section will discuss the overall purpose and 
general operation of preprocessors and processors in the ESi 
system. A more complete discussion will be found in Reference 
20 . 

a. Preprocessor 

Basically, the preprocessor must prepare the 
signal received from the receiving system for preliminary 
analysis and for advanced analysis in the main processor . 
Figure 2.8 is an ESM preprocessor. [20; p.83 

After digitization, the initial parameters of 
pulse width, time of arrival, direction of arrival , amplitude, 
and frequency must be measured on a pulse by pulse basis. 

<See Parameter Measurements in this chapter.) 

The pulse data must be compared with previously 
received and analyzed signals to see if this signal has 
already been processed. Pulse data is compared with signal 
data in the files that contain the active emitters and the 
uninteresting emitters. If a match is present, then the data 
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Figure 2.8. Preprocessor 
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can be deleted from the process. Reduction o-f the data stream 
and sorting o-f the signals is the purpose and result o-f this 
deletion in addition to avoiding unnecessary processing. One 
o-f the main Jobs o-f the preprocessors is to convert the 
approximately 10^ to 10® pulses per second being received to 
about 1000 pulses per second that a good main processor can 
handle. C21: p.1^43 

Unmatched or new signals must be dei n ter 1 eaved 
(see deinterleaving in this chapter) and processed in the PRI 
processor section. Calculation o-f the PRI is usually done 
simply by subtracting TOA''s o-f similar pulses. Staggered 
PRI's and jittered PRI's have to be handled with a more 
complex algorithm. [20: p. 121 

Ni th a known PRI, TOA prediction, additional 
sorting, and comparison o-f calculated PRI's is done until a 
PRI is con-firmed to exist in the data. This in-formation, the 
newly calculated PRI and the initial parameters o-f the pulse, 
is then passed to main memory o-f the processor -for analysis 
and updating o-f the active -files so that pulses o-f the 
incoming stream can be deleted as soon as possible. 

In summary, the preprocessor prepares the incoming 
signals -for processing, measures their initial parameters, 
per-forms preliminary analysis and sorting, and passes the data 
in-formation to the main processor. 



47 




Figure 2.9. 



Main Processor 
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b. Main Processor C20; pp. 13-133 

A-fter all the action in the preprocessor and the 
data reduction that has been accomplished, there is still 
plenty to do. The -functions o-f the main processor are: 

1. Advanced analysis o-f the signal for scan rate and exotic 
emitter <such as chirp radars) identification 

2. Emitter and platform identification 

3. Updating the emitter files 

4. System control and monitoring 

5. Interfacing with the operator and ECM equipment 

Since the processor is usually a general purpose 
computer, additional functions of navigation or routine system 
calculations could be added. It can also help the PRI 
processor handle its calculations. 

Emitter identification is accomplished by 
comparing the measured emitter parameters with emitter 
parameter lists until an identification is made or the signal 
is tagged as unidentified. Interface with the operator for 
manual analysis is often the result of an unidentified 
signal . 

3. Display System 

An operator must have a visual presentation of the 
signal and its parameters. Almost universally CRT displays 
are used. CRT's usually have an intensity modulated or a 
deflection modulated design, with advantages and disadvantages 
in each one. 

Multi-gun CRT displays can create a combined display 
that offers a pulse analyzer, presenting frequency, PRF, PW, 
and the presentation of an amplitude spectrum. 
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A warning system usually operates in conjuction with 
detection o-f high threat signals. Flashing lights and/or an 
aural alarm will occur, indicating immediate attention is 
required. 

The audio system operates in parallel with the warning 
system. 0-f-fering not only alarms to high threat signal 
reception, the system usually has accomodations that allow the 
operator to monitor the received signal through headsets or a 
speaker. A trained and experienced operator can determine 
scan rate and analyze the scan modulation type i-f the signal 
density is low enough. 

Bearings -from DF equipment are usually presented on a 
polar display. The direction is usually indicated in relative 
bearings o-f 0 to 360* around the plat-form, with some output or 
means o-f determining the true bearing -for plotting purposes. 

Computer management o-f an ESM system allows for 
most calculations of bearing, PRF , PW, and frequency to be 
accomplished and shown as soon as the bearing is displayed to 
the operator. An excellent example of ccxnputer managed 
displays is the ALQ-78 ESM system on the U.S. Navy P3-C 
Orion. Digital readouts of the PRF, PW, and frequency are 
available next to the displayed bearing from the aircraft. 

The bearing remains on the screen and updates automatically 
until removed by the operator or the system. With several 
cross bearings, the emitter location can be obtained since the 
computer is not only managing the display but the navigation 
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o-f the aircraft. By "hooking" the place where the bearings 
cross, the operator can cause the computer to calculate the 
1 at/1 ong of the emitter fix. 

B. ESM SIGNAL PROCESSING 

The acquisition of signals is usually a design problem 
that is driven by the character isti cs of the signals to be 
acquired. The designer chooses antennas and receivers to 
detect signals with regard to the mission, the platform, and 
the frequency band of interest. Generally, he desires the ESM 
system should have rapid acquisition, high sensitivity and 
frequency resolution, and wide frequency coverage. 

After these signals are acquired, signal processing is the 
order of the day. Processing is required immediately to 
determine the high-threat signals <especial1y for tactical 
aircraft) , and a lot of processing is required in today''s 
dense electromagnetic environments. The signals are getting 
more complex also, with jittered and staggered PRF's and 
variable PW's in addition to frequency agility and new 
schemes of radiating and detecting a return. 

Computational speed becomes one answer to solve these 
processing riddles, and the debate between digital and analog 
techniques fuels upon this quest for speed. 

Surface Acoustic Wave <SAW) devices and Charge Coupled 
devices <CCD's) are being challenged by digital technology 
staked by the Pentagon's Very High Speed Integrated Circuits 
Program <VHSIC) . These digital IC's are being structured to 
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per-fopm high speed calculation, with a N X N matrix inversion 
(a N3 operation) being a gauge o-f successful application. 
Today's analog technologies currently offer high speed also, 
but will require more precision and dynamic range to meet this 
challenge. L22: p.?l] 

In the end, the ESM system designer will win as new 
digital arch i tec tures and new device structures will yield 
increase in speed and throughput. 

This section will discuss important signal processing 
parameters of the ESi system, and the techniques used to 
separate the signals for identification and processing from 
the immense amounts of pulses being received. The entire 
process is a continuous flow of calculations, comparisons, 
and decisions involving tolerances in frequencies, pulse 
width's (PW> , pulse repetition intervals (PRI) , and other 
parameters . 

1 . Parameter Measurement 

Listed below are the initial parameters, as provided 
by the receiving system to the ■ preprocessor . They are the 
basic tools that are used to deinterleave < see next section) 
the pulses and match them to existing tracked or not tracked 
emi tters. 

1. Direction of Arrival 

2. Ampl i tude 

3. Frequency 

4. Pul se Width 

5. Time of Arrival 
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Those that aren't matched must be analyzed further by 
calculating a pulse PRI and/or undergoing advanced analysis. 
Calculation of the PRI is done in the preprocessor . 

Advanced analysis is accomplished in the main processor using 
derived higher order parameters during the task. 

A discussion of scxne aspects of how the initial 
parameters and higher order parameters are measured follows in 
this section. [21: p. 1663 

a. Direction of Arrival <D(^) 

An important parameter, the DOA could be measured 
with a rotating DF antenna but this greatly lowers the 
Probability of Intercept <POI) . A more usual method uses a 
multi-port (antennas) bearing measurement subsystem that has a 
high <%100%) POI . 

Phase comparison techniques can also be used, 
providing greater accuracy but with a tradeoff of smaller 
bandwidth . 

b. Ampl i tude 

Usually the amplitude is termed by the peak value 
of the received pulse with some tolerances included to provide 
for variances in emitter radiated power. 

The amplitude parameter is used in the advanced 
analysis section of the typical processor to determine a scan 
pattern, but is not considered a good sorting parameter 
because of propagation corruption, multi-path, and other 
variables that influence the amplitude of the received pulse. 
C16: p. 2g-33 
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c. Frequency 

Frequency may be measured with a scanning 
superheterodyne receiver, but this method has a low POI -for 
reasons similar to those o-f using a rotating antenna -for DF . 

A channelized receiver provides a frequency measurement but 
this usually is a much more expensive system. 

Most common is the IFM receiver (see section on 
receivers in this chapter), which has an excellent resolution 
and high POI . 

d. Pulse Width 

Pulse widths can be measured by simply noting 
when the pulse rises and falls through a threshold. 

Reflections tend to severely corrupt this method of 
measurement, however. 

A more accurate method is to sample the pulse many 
times and take the PW in the distance between the 3 dB 
points. 

e. Time of Arrival <TOA) 

Threshold crossing measurement of this parameter 
could also be used as in the case of the PW, but a better 
measurement can be obtained by noting the TOA of the first 3 
dB point. TOA's are used in the calculation of PRI's. 

Adjacent pulse TOA's are simply subtracted to find the PRI . 

f. Higher Order Parameters 

Parameters that are calculated from the other 
initial parameters are the derived or higher order parameters. 
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The PRI is used to link interrupted or spl i.t 



chains o-f pulses -from the same emitter. It can also be used 
as a parameter to identify a jittered or staggered PRI. 

Scan patterns can be recognized by analysis o-f the 
pulse amplitudes generated as the radar beam sweeps across the 
plat-form. Knowledge o-f the scan pattern can be o-f help in 
identi-fying the emitter since the radiation pattern can be 
correlated to emitter type. C21: p. 1671 
2. Deinterleaving 

A stream o-f pulses -from several di-f-ferent emitters 
includes di-f-ferent frequencies, amplitudes, and PRI's mixed 
together (interleaved). Separating these different pulses 
into chains of similar pulses for identification is the 
deinterleaving problem. 

With high density environments, the time available for 
processing each pulse may be only around 1 microsecond. C21: 
p. 168] Therefore, deinterleaving algorithms, if done 
digitally, must be fast. Analog equipment that sort pulses 
have the same requirement. This required speed is the 
obstacle that hinders the use of Digital FFT's for pulse 
sorting. 

Provision must be made in algorithms for missing 
pulses in the data stream. Missing pulses is one of the 
greatest confusion factors in the system. Methods are 
currently in use to smooth fluctuations of the number of 
pulses that are input into the preprocessor . C21; p.l67] 
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General discussions o-f deinterleaving methods are 
given in this section. The -first two are digital methods, 
using the measured parameter o-f the pulse, and use o-f one or 
the other depends upon the accuracy of the E3i system. Two 
hardware methods are briefly covered also. 

a. Cell or "Pigeon Hole" Techniques 

A system that accurately measures the parameters 
of the pulse stream can use this technique as a fast and 
efficient deinterleaving scheme. 

Pulse data with similar parameters are directed 
into cells or "pigeon holes". The cells soon contain alike 
parameters, such as bearing, frequency, or pulse widths that 
can be analyzed for emitter identification or further sorting 
if too wide of a variation exists. 

Variations could be caused by emitter variations 
itself or purposeful variations such as a jittered, staggered 
PRI , or frequency agile radar. If variations (unmatched 
parameters) are present, then they are passed to a more 
complex algorithm for processing for the special modulations. 
Otherwise, a simple analysis or comparison with the active 
emitter file or emitter parameter file results in an emitter 
identification . 

b. Time of Arrival Techniques C21; p. 168] 

Nhen accurate computation of pulse parameters 
isn't available, a time slice of the pulse stream is used to 
simplify the data under analysis. A time slice will 
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contain about 8 to 512 pulses. This scheme reduces the number 
o-f pulses being considered in a high density enviroment, and 
■forces comparison o-f parameters to those only in the present 
time si ice. 

The -first pulse in the time slice is used as a 
re-ference and the rest o-f the pulses is scanned for a match. 
When a match is found, the TO^ of the pulses is used to 
calculate a PRI . When a sufficient number (usually 6 or 7) of 
matching pulses are found, then the data is grouped as a chain 
and is output to the main processor. A problem generated with 
this technique is the splitting of chains in the pulse data 
stream by the time slice. Several time slices may have to be 
analyzed before enough matches in a slice are found to form a 
chai n . 

c. Analog Methods 

A pulse sorter operating on the pulse period 
that has some analogy to digital TOA methods is currently in 
use. C 16: 2g-7 to 2g-9] Consider Figure 2.10. 

The pulse stream is fed into a coincidence 
detector-discriminator that attempts to match the pulses with 
gates generated by a gating pulse generator. The gating 
pulses are initially generated at a PRI slightly higher than 
the highest pulse PRI. By gradually reducing the gate PRI, 
the highest PRI pulse train will be captured by the gates if 
the phase difference between the gates and pulses of the 
stream. This phase difference is used to generate an error 
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signal that causes the gate PRI to decrease until synchronized 
with the highest PRI pulse train o-f the stream. The pulse 
train is n«g sorted out of the stream, and the remaining 
pulses are applied to a second coincidence 
detector-discriminator where the process is repeated. 

Other methods are used on the pulses left in 
the stream after passing throught all the 

detec tor— di scr imi nators . Staggered PRI"s can also be handled 
by using two or more gate generators operating in 
synchronization . 

d. Frequency Dariain Methods 

Bandpass comb filters C16: p.2g-91 can be 
used to pass only harmonically related frequency components. 

By tuning the filter with the pulse stream as an input, the 
filter will pass components when it is tuned to a PRI of the 
pulse stream. Disadvantages of this method include the low 
resolution of the sorted PRI's and inability to distinguish 
harmonic components of different PRI's. 

Digital FFT's are attractive as a sorting 
tool, but hardware (size and cost) and computational speed 
have a way to go before it offers practical alternatives. 

VHSIC will certainly improve both the hardware and speed 
aspec ts . 
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III. THEORY OF NALSH AND RADEMACHER FUNCTIONS 



Electrical Engineering could be said to revolve around the 
sine and cosine -functions. They are the basis -for development 
in many areas due in part to the inherent properties in their 
■frequency domain representations. 

The application o-f digital techniques and semiconductor 
technology, though, has brought -forth uses and awareness o-f 
other orthogonal -functions. These o-ften do not have the 
desirable properties o-f sine and cosine -functions -for use in 
linear, time invariant networks, but they do have other 
advantages that render them use-ful in other applications. 

The early part o-f this century saw the introduction of 
several two level orthogonal functions, or sometimes called 
binary functions because of their amplitudes taking only two 
values. Work in orthogonal matrices, by Sylvester in 1367 
and Hadamard in 1893, was an early approach to these 
functions. 

In 1910, Alfred Haar presented a complete set of 
rectangular functions that took only two values but yet 
provided complete expansion of a continuous function. These 
Haar functions could be made to converge uniformly and 
rapidly. See Reference 23 for additional information. 

A German mathematician, H. Rademacher , presented another 
set of orthogonal functions in 1922. These were followed by 
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the Walsh Functions, defined in 1923 by the American 
mathematician J.L. Walsh. The Rademacher -functions, although 
independently presented, were -found to be an incomplete but 
true subset to the Walsh -functions. C23; p.v-vi] 

This chapter will brie-fly examine the Rademacher 
•functions. They are important because products of the 
functions yield a certain ordering of the Walsh functions. 

The Walsh functions will then be examined in more detail. 

The reader is referred to Appendix B for a generic treatment 
of orthogonality and orthogonal functions. 

A. RADEMACHER FUNCTIONS 

These two-level orthogonal functions are represented by 

R<n,t) <3.1) 

n = 0,1, 2, 3,... 

and can be seen in figure 3.1 for n = 0 to <6. C23: p.6I 

Notice that they are a series of rectangular pulses or 
square waves and have 2*^“^ periods over a time base from 0 
to T. Mathematically they could be defined by 

R(n,t) =sign<sin<2‘^ltt)> <3. 2) 

The first Rademacher function, R<0,t>, is equal to one for 
the entire interval, and subsequent Rademacher functions have 
odd symmetry with amplitudes of -*• 1 and -1. All have unit 
mark-space ratio. 

One can generate these functions with a sine function of 
appropriate frequency with amplification and hard limiting. 



61 



The appropriate frequency would be one with the same zero 
crossing positions as the Rademacher functions. C23: p.7] 




R(0,t) 



Rd.t) 



RI2.I) 



R(3.i) 



R(4,t) 



R(5.l) 



0 T 

Figure 3.1. The First 6 Rademacher 



Functions. 



B. WALSH FUNCTIONS 

Of more importance than the other two-level orthogonal 
functions discussed are the Walsh functions. They form a set 
of rectangular waveforms taking only two amplitudes, +l and 
-1, but do not have a unit mark— space ratio like the 
Rademacher functions. They are defined over a time interval 
T, which must be known to assign values to the functions. 

Two arguments are required for complete specification of 
a Walsh function 

WAL<n,t) (3.3) 
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The time base Ct] of the function is usually specified as 
t/T, and thus normalized from 0 to 1 . The number Cnl is 
equal to the number of zero crossings a Walsh function has 
during the time base. C23: p.7] 

Figure 3.2 shows the Walsh functions. The syiwnetry of 
the functions and the concept of sequency (next section) was 
used by Harmuth C24] to define another notation for each 
Walsh(n,t). Each WAL(n,t) is either odd or even about the 
midpoint. If the function is odd, then it can be referred to 
as SAL(k,t) . If the function is even, then l«WL(n,t) is also a 
CAL(k,t). Note the similarities between this notation and the 
sine-cosine functions. Relationships between Cnl and Ck3 are 

CAL(k,t) = WAL(2k,t) (3.4) 

SAL(k,t) = WAL(2k-l,t) 

This notation, SAL(k,t) and CAL(k,t) , is seen on the right 
side of Figure 3.2. 

1. Sequency, Ordering, and Phasing 

Sequency is a term proposed by Harmuth to describe a 
periodic repetition rate which is independent of waveform. It 
is defined as "one half of the average number of zero 
crossings per unit time interval* and abbreviated "Zps". C23: 
p.l31 A 100 HZ sine wave has 200 zero crossings per second, 
so frequency is a special measure of sequency as applied to a 
sinusoidal waveform. 

Ordering the functions by ascending number of zero 
crossings (WAL(n,t) notation) is called sequency ordering. 
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Figure 3.2. The Mai sh Functions in Sequency OrcJer. 
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Although Cn] is not the sequency o-f the -function, the 
equivalence o-f the NAL<n,t) with SAL<k,t) and CAL<k,t) 
notations automatically arranges the -functions in ascending 
order o-f sequency. The value o-f Ck] in the SAL and CAL 
function is the sequency for the function. Sequency ordering 
is the preferred ordering for spectral analysis and filtering. 

The natural ordering (or Pal ey ordering) is obtained 
by generating the Walsh functions with products of Rademacher 
functions. Figure 3.3 shows the Walsh functions in natural 
order, which are referred to as 

PAL(n,t) (3.5) 

n — 0,1, 2, 3,... 

because Pal ey C243 used this ordering. For theoretical and 
mathematical work, image transmission, and computational 
efficiency, this ordering is usually the preferred one. C23; 

p. 18] 

A modulo— 2 addition (Appendix C) relationship exists 
between the Walsh functions with sequency ordering and the 
functions with natural ordering. C23; pp. 31-32] An example 
will illustrate the conversion between the two orderings. 
Consider WAL(n^g,t) = WAL(13,t). 

Here n^^ = 13 = 11012* To find the equivalent 
natural order function, PAL(np,t), first add a zero (base 
2) as the leftmost digit to base 2 n^,^. Then generate the 
value of Op in base 2 by doing modulo— 2 addition between 
consecutive zeroes and ones of base 2 n^j, starting with the 
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Figure 3.3. The Walsh Functions in Natural Order 
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le-ftmost digit <the zero that was added on to base 2 n,^> . 
Following these rules, = 13 = II 0 I 2 = 0110l2> and 

0 ® 1 = 1 = K3 

1 ® 1 = 0 = K 2 

1^0= 1 = kj 

091 = 1 = k0 

So , 

np — (k3k2l^ 1^0^2 ~ 101l2~ 

Figure 3.3 shows the result is correct. 

WAL< 13, t) = PAL< 1 1 , t) 

The phasing o-f the functions as given in Figure 3.2 
and 3.3 is known as Harmuth C25] phasing. This phasing 
emphasizes the phase similarities with sine-cosine functions. 

t 

Positive phasing C23: p.l8] is where all of the 
functions begin at +1. This involves a sign change for many 
of the Figure 3.2 and 3.3 functions, and is the result of 
function derivation from Hadamard matrices (orthogonal 
matrices composed only of + 1 and -1 elements; see Ref. 23, 
pages 24-25) or Rademacher functions (next section). 

2. Derivation of the Hal sh Functions 

Products of selected Rademacher functions will yield a 
complete set of Hal sh functions in natural order. Recall that 
the natural ordering of the Hal sh functions was referred to as 
PAL( n , t) . 

In terms of Rademacher functions, 

m+ 1 

PAL(n,t) ='"JP b^R(i ,t) (3.6) 

i=l 

where bj is either a zero or a one, indicating the presence 
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or not o+ the i ' th Rademacher -function. To -find the bj's 
and Cm], represent Cn3 as a binary number, 

with m = the highest power o-f two -found in the binary number 

Cn]. The bj's present in the positions o-f base 2 Cn] 

indicate the presence or not o-f the Rademacher -function in 

Equat i on (3.6). 

Consider PAL(9,t). Then 

n = 9 
n = 100 I 2 

with m = 3, b^ = 1, bj = 1, and -finally, 

PAL(9,t) = R(4,t)R< l,t) 

It should be noted that this product is positive 
phasing. Harmuth phasing can be obtained by de-fining the 
Rademacher -functions over the interval -J4<t<^. This 
de-finition inverts all o-f the -functions shown in Figure 3.1. 
Products o-f the inverted -functions will yield natural ordered 
and Harmuth phased Mai sh -functions. C23: p.223 The Walsh 
-functions with sequency ordering can be derived with 
Rademacher -functions by use o-f the Gray Code (Appendix C) . 
Consider WAL(13,t) where n = 13. 

In Gray Code, n = 13 is represented by 1011, and 
recognizing the ones o-f the Gray Code are in the -fourth, 
second, and -first positions, (similar to the derivation 
o-f the PAL(n,t) -functions), C23; p.23] 

WAL(13,t) = R(4,t)R(2,t)R( l,t) 
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other derivations o-f the Walsh -functions can be 



accomplished by di-f-ference equations, Hadamard matrices, and 
Boolean synthesis. See Re-ference 23, pages 28-26. 

3. Walsh Series 

It is well known that a time -function f<t) can be 
expressed as a sum o-f a series o-f sine and cosine -functions. 
Each -function has a coe-f-f icient that determines the value o-f 
the function in that series. 

The time function can be expressed in a similar way 
using the Walsh functions. C23: p. 131 



N-1 

f(t) = a 0 WAL( 0 ,t) + Z a„WAL<n,t) 

n=l 

where the coefficients would be calculated frcwn 



<3. 3) 



a0/2 = 1/T 



i: 



<t)WAL<0,t)dt 



an = 1/T 



y 



< t)WAL<n , t) dt 



(3.9) 



<3. 10) 



with N = the desired number of Walsh terms, and T is a given 
length of the function. 

Each Walsh function has a coefficient associated with 
it that gives the value of the function in the series the same 
way as does the Fourier Series coefficient. [23: p.l3] 

If the function is periodic, and T is normalized to be 
equal to 1 , then a 0/2 is the mean of the function, since 
WAL<0,t) is equal to 1 over the interval 0 <t<l. The 
expressions for a 0/2 in both the Walsh and Fourier Series 
are equal with this normalization. 
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Consider that the -function is absolutely integrable 
in the interval 9<t<l. Then the -function can be expanded in a 
Walsh series o-f the -form [26: p.232 and 23: p.40] 

f<t) = Z a<n>WAL<n,t) (3.11) 

n =0 



where 



a< n) 



■f < t)WAL(n , t) dt 

0 



<3. 12) 



A more general de-finition -for a periodic -function 
an interval <0,T) is 

04 



a< n) 



•f(t) = Z a(n)WAL(n,t/T) 
n =0 



= 1/T 




t)WAL(n,t/T)dt 



de-fined over 



<3. 11a) 



<3. 12a) 



Now coe-f-f i ci en t a 0 becomes the mean o-f the -function, 
since WAL(0,t/T) = 1 over the normalized interval 0St/T<l. 

The Walsh Series can be de-fined over the same interval 
<0,T) using the CAL and SAL notations, and an in-finite number 
o-f terms: C26: p.2321 

w w 



-f(t) = a0WAL(0,t) -*• Z 

i=l 



Z (aiSAL(i,t) + bjCAL(j,t)) 
j=l 

<3. 13) 



wi th 



a 5 <n)= 1/T 



a^.(n) = 1/T 



$:■ 

S; 



<t)SAL(n,t/T)dt 



<t)CAL(n,t/T)dt 



<3. 14) 



(3.15) 
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A couple o-f simple examples follows to illustrate the 
Walsh Series expansion. 

a. Expansion of Sin<t) in a Walsh Series 

The sine function will be defined over one period, 
with a period of 2'K, and the Walsh functions will be defined 
as t/T goes from 0 to 1 . 

Calculating the first coefficient, using Equation 

(3. 12a) , 



a 0 = 1/21C • \sin< t)WAL<0 , t/T) dt 
J 0 



^0 

^0 



p2K 

1 / 2 K • \sin(t) ( Ddt 
/0 



= 0 



This is a confident calculation, 
Know that the average of a sinusoid is 0‘ 

The second coefficient, 

p2lC 

a^ = 1/21C • \ sin(t)WAL(t,t/T)dt 

J 0 



pi 



ai = 1 / 2 H • C\sin<t) <-l)dt + 

J0 

aj = -2/X = -.63611 



5 



A lot of menial integration will 
coefficients of any practical value; 



as most students 



2lt 

sin<tXl)dt 3 
1 C 



yield only these 



as = .26348 
a ^2 “ • 12653 

Coefficients numbered 2,6,9, and 10 are present 
but have values less than 0.06. 
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The -function sin<t> is then approximated by 
sin<t> = aiUJAL<l,t) + agWALCSjt) + ai3WAL<13,t) 
where the coe-f-f icient values are given above. 

Figure 3.4 C23; p. 15] shows the three Walsh 
■functions used and their sum. The negative sign o-f the 
coe-f-ficient aj has been used to invert the Walsh<l,t) 
■function . 





Figure 3.4. Sin<t) with a Limited Walsh Series 
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Of course, the addition of more of the terms of 



the Series would yield a better approximation. One can see 
that representing a smooth curve using rectangular shapes 
would require a large number of terms of the Mai sh Series. 

The Walsh Series will give better results 
representing a rectangular function. 

b. Expansion of a Rectangular Function 

A rectangular pulse function should lend itself 
nicely to representation by other rectangular functions. 

Indeed, it is hoped that the representation might 
be exact or with small error in only a few terms. 

Consider this function; 



f<t) = 1 



0<t<.25 

25<t<.75 

75<tn.0 



0 

1 



Proceeding with the coefficients 




f(t)WAL(0,t)dt 

0 




0.5 



and 



a2 



25 

<1X-I)dt + 
0 



I 



1.0 

< 1) <-l)dt 
75 



5 



Coefficients a^, ag, and above will all 



calculate to be zero 
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This rectangular time -function, uahich could 
represent two radar pulses, can be represented exactly by two 
Walsh -functions. 

■f<t) = 0.5WAL<0,t) - 0.5WAL<2,t> 

Figure 3.5 shows the -function, the Walsh 
-functions, and the sum o-f the Walsh -functions. This is quite a 
di-fference -from the representation by Fourier Series, which 
would require a lot o-f terms to present an inexact 
represen tat i on . 

c. Wave-form Synthesis 

These results obtained above indicate that a 
continuous waveform is more suited for Fourier transformation 
and a discontinuous waveform, certainly a rectangular one, is 
more suited for Walsh transformation. This conclusion should 
be easy to accept with examination of the structure of the 
orthogonal functions <sine-cosine and Walsh). 

Figure 3.6 C23: p. 333 shows a reconstruction of 
two waveforms and the number of Fourier and Walsh terms used. 
The reconstruction supports the conclusion very well. In 
addition, Beauchamp [23: p.353 reports that the same 
conclusion can be drawn when the Fourier and Walsh Transform 
are used to reconstruct a continuous and a step waveform. 

In coming to this conclusion, Beauchamp 
considered the effect of the number of allowed levels for 
quantizing the waveform and the sampling interval. C23: p.353 
The sampling interval determines the number of coefficients 
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Figure 3.5. Expansion o-f a Rectangular Function 
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Figure 3.6. Wave-form Synthesis Using Walsh Functions. 
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produced, and the highest Walsh coe-f-f icient found from the 
Transform determines the number of data points available for 
any quantized level. This is a built-in limitation on the 
accuracy of the Walsh Transform that the Fourier Transform 
doesn't have. 

d. Digital Sampling 

A theorem similar to the Sampling Theorem for 
sine-cosine functions is applicable to the Walsh functions. 

The minimum sampling rate is C23; p.373 

fs = 2 ^'*^ <3. 16) 

with k being the power of 2 that represents the sequency 
bandwidth. The sequency bandwidth would be the sequency of 
the highest component desired to be represented. 

Suppose a CAL(12,t) ccnnponent had been determined 
and it was desired to include this component in the bandwidth. 
The sequency of this component is 12, and the closest power of 
2 that is at least 12 is 4 <2^ = 16) . Thus k = 4 and 

fs = 2^+1 = 32 

For canparison, a cosine of frequency 12 would 
require a sampling rate of 24 <23Sfrequency) . Aliasing is 
equally applicable to the Walsh functions and should be 
considered. C23: p.373 

4 . The Walsh Transform 

Section 3 stated that a continuous function over the 
interval 9<t<l could be represented by a sum of Walsh 
functions, shown in Equation 3.11. 
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■f <t) 



<3.11) 



CO 



= Z ar,l4i»L<n , t) 
n=0 



where 



- ■ i 



< t)HAL<n , t) dt <3.12) 

I 

Thus, a trans-form pair can be de-fined: C23: p.40] 



w 



■f < t) = Z F<k)WAL<k , t) 
n=0 



<3.17) 



K> =j‘ 

J0 



F<k) =\ -f<t)WAL<k,t)dt <3.18) 

/0 

F<k) could be written, -for an interval 0<t<T, 



f. 



F<k) = 1/T ] -f < t)WAL<k , t) dt 



<3.19) 



Now suppose the interval T is divided into 
N parts, letting t = i^, where ^ is a small time increment and 
i is an integer. The aim is to convert the integral into a 
summation -for discrete computation. 

Then , 



f 



wi th 



F<k) = 1/T \ -f<t)WAL<k,t)dt 



N-1 

S5 1/N-d • Z -f < i.d)WAL<k , iiJ) 
i=0 



-d = dt 
= T 
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Remembering ^ to be a given small increment, but 



unchanging in value, 

N-1 

F<IO = 1/N • Z -f<i)WAL<k,i) (3.22) 

i=0 

wi th 



•f(i) = evaluation of •f(t) at intervals iA 
— a series o-f numbers 
Change the notation -from k to n, -f<i) to 
xj , and F(k) = , and write the -finite discrete Walsh 

trans-form pair. 

N-1 

X„ = 1/N • Z XjWAL(n,i) <3.23) 

n=0 



n = 0 , 1 ,2, . . .N-1 

and, 

N-1 

xj = ZX„WAL<n,i) <3.24) 

n=0 



i = 0 , 1,2, . . .N-1 

The CAL and SAL trans-forms -for a discrete series can 
also be -found by using the equivalence of WAL with CAL and SAL 
notations; 



N-1 

Xc<k) = 1/N . Z XjCAL<k,i) <3.25) 

i=0 

N-1 

X^.<k) = 1/N • Z XjSAL<n,i) <3.26) 

i=0 

An even or odd sequence when trans-formed will have 
similar properties as the discrete Fourier Coe-f -f i c i en ts -from 
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even or odd -functions. That is, an even sequence (symmetrical 
around its midpoint) will trans-form into only CAL -function 
coe-f-ficients, and an odd sequence will trans-form into only SAL 
coe-f-f icients. C23: p.413 

The Discrete Walsh Trans-form (DWT) has some advantages 
over the Discrete Fourier Transform (DFT) because it involves 
only additions and subtractions, not multiplications as in the 
DFT. Also, the transform is not noisy because precise 
representations in the digital computer is possible. Precise 
representations of sine-cosine functions in the DFT is not. 
Finally, the Walsh transform kernel is ±1 and is considerably 
easier to calculate. C23: p.423 

An important point of the transform is that the CWT is 
its own inverse, and a separate inverse transform definition 
is not required. See Equations (3.23) and (3.24). 
a. The Fast Walsh Transform 

Computation of the DWT would involve N2 
mathematical operations of either addition or subtraction. 

This number of operations has been decreased by using the 
redundancy in the Walsh matrix representation of the Transform 
to N1 og 2 N operations, the same number as in the DFT. C23; 
p.543 Remember that the DWT involves additions/subtractions, 
not complex multiplications/additions as in the DFT. 

The Fast Walsh Transform can be described in a 
signal flow diagram that results in a “butterfly" appearance 
just as in the FFT . See Figure 3.7 on the next page. Unique 
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Figure 3.7. FWT "Butter-fly"; N = 16. 
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■features o-f the butterfly in the final stages of calculations 
offer an opportunity to calculate the FWT "in-place", without 
using computer working space. C23: p.55] The result is not in 
sequency order, however, and a sorting routine must be 
included if sequency order is desired. 

Appendix A gives listings of the FWT programs used 
in this thesis. 

b. Walsh Coefficients and the Walsh Matrix 

The computation of the coefficients from the FWT 
is very fast and takes advantage of the properties of a 
Hadamard matrix C23: p.52] 

The Walsh matrix is related to the Hadamard matrix 
and can be used to determine the coefficients of the Walsh 
Series. C27; p.4] Consider this matrix equation, 

Ca^] = 1/N • CW] . Cxi] <3.27) 

where tWl is a N X N matrix of the Walsh functions, Cxj] is 
a column vector of N sampled values of the function, and 
Cap], n = 0 to N-1, is a column vector of the Walsh 
coefficients. 

The Walsh matrix can be thought of as the matrix 
of sampled Walsh functions, and thus, is composed of plus and 
minus ones. To compute the coefficients, simply perform 
addition and subtraction of the sampled values according to 
the row of the Walsh matrix being used. An example will show 
this computation. 
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Let the number of samples in the column vector = 3 



= N. Nriting the matrix equation, omitting the 1/N term, 

<3. 16) 



= tW] 



CXi 3 



^0 




^1 




a2 




^3 




a4 




^5 








^7 





+ + + + + + + + 

+ + + + -- -- 
+ + -- -- + + 

+ + -- + + -- 

+ -- + + -- + 

+ -- +- + + - 

+ -+ -- + - + 

+ - + - + -+ - 



For comparison, use N sampled values of the 
rectangular function on page 73. Thus, 

=[ 11000011 ] 

i = 0jlj2j3j««»j7 





5<0 




^1 




^2 








X4 




^5 




X6 




[x7j 



Now evaluate the constants a^ using matrix 
multiplication. A plus sign in the Walsh matrix equals +1, 
and a minus sign means -1. Don't -forget to divide each 
constant by N. 

a0 0.5 0 

aj = 0 a5 = 0 

a2 = 0 • 5 a^ = 0 

a0 = 0 a7 = 0 

Comparing with the value o-f a 2 given in the 
Walsh Series expansion o-f this same -function (page 73> , we 
note a sign change. The reason lies in the use o-f positive 
phasing -for the Walsh matrix and the use of Harmuth phasing in 
the Series expansion computation. 

The same answer has been reached, although this is 
the discrete case. 



83 



I 



Matrices can also be used to show the Inverse 
Discrete Walsh Trans-form. From Equation <3.24), 

N*[Xj] = EW] -Ca^] <3.23) 



= CW] -CW] -Cxj] 



and 



[Xj] = 1/N-CWI -CW] -CXil <3. 2?) 

Equation <3.29) says to get the sampled values o-f 
the original within a constant N, just multiply the 
coe-f-f icients by the Walsh matrix CW] . 

c. Input Time Shi-fts and Their E-f-fect 

The DFT magnitude is invariant to the phase o-f the 
input signal. Although Walsh Trans-form signals con-form to 
Parseval's Theorem, the CMT is not invariant and the same 
spectral representation cannot be achieved independently o-f 
the phase or time shi-ft o-f the signal . E23; p.42] 

The DWT is time invariant when the time shi-ft is 
obtained by dyadic translation C23: p.433 . Namely, i-f x^, 
the input number series, is shi-fted to Zj = Xj^, where 

i®p = modulo -2 addition 

then , 

2c^<K) + 2^2<k) = X(-2<K) + X^2<K) 

i-f the coe-f-f icients are expressed in CAL notation <Z,-, X,-) 
and SAL notation < 23 ^, X 3 ) . See Appendix C -for modulo-2 
addi tion . 

The importance o-f variations with time shi-ft can 
be lessened by considering sums o-f squares o-f coe-f-f icients. 
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This is one de-finition o-f the Mai sh Power Spectrum, and is 
considered in the next section. 

5. The Walsh Pcx«jer Spectrum 

Spectral analysis using sequency has some 
advantages over analysis in the -frequency (Fourier) domain. 

A de-fined pcxaer spectrum -from the Walsh -functions is easily 
applied to discontinuous or time limited functions, for 
example, and it is possible to have a sequency limited 
spectrum for a time limited function in contrast to the 
Fourier Spectrum. C23: p.871 

Although there are several derivations of the 
Walsh pcwer spectrum C23: p.?l], a more normal, i .e. 
analogous to Fourier, approach will be used to make a spectrum 
definition. As there are readily available methods (R4T) for 
calculating the Walsh coefficients, a combination, 
specifically the sum of the squares of related coefficients, 
would give an easy to use definition. Also, as stated 
earlier, by using the sum of the squares of certain 
coefficients, the time shift variance that plagues the Walsh 
Transform can be lessened. 

Define, then, the Power Spectral Density (PSD) 
coefficients of the Walsh spectrum. C23: p.l083 

P(0) = X^.2(0,t) (3.31) 

P(k) = X^2(k) -»• X52 <k) 

P(N/2) = X 32 (n/ 2 ) 

k = 1,2, .. . ,(N/2 - 1) 
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Note that (N/2 + 1) spectral points are 
generated. This method is analogous to Fourier power spectral 
analysis, where the power coe-f -f i c i en t is the sum o-f the 
squares o-f the real and imaginary parts o-f the complex Fourier 
Trans-form coe-f-f icients. 

Beauchamp C23; p.l00] states that taking the 
square root o-f P<k) does not give a sequency amplitude 
spectrum, but other re-ferences C26: p.233] use -4P|^ as the 
de-finition o-f the sequency spectrum. 

Another spectrum can be de-fined using the Walsh 
Coe-f-f icients. The name Group Spectrum is the author's term 
and will be used here. Beauchamp [23; p.1061 re-fers to it as 
"the odd-harmonic sequency spectrum", and Campanel 1 a/Robinson 
C26: p.234] re-fer to this spectrum de-finition as a generalized 
Wal sh-Four i er spectrum [283. 

The power content o-f individual groups o-f 
-frequency components has been shown to be equal to groups o-f 
components in the sequency domain. Any group o-f components, 
in either domain, is a collection o-f the orthogonal components 
that make up the original signal or -function. 

Call groups o-f sequency components G^ , with [n] 
an integer, and consider a discrete signal made up o-f a 
sequence o-f length N. 

The total energy o-f the signal is, by Parseval's 
Theorem, equal to the sum o-f the orthogonal components [293 . 
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<These orthogonal sequency components have been grouped into 



6 ^.) Then, 

N- 1 m 

1/N • 2sy 2 = ZGn 
k =0 ' n =0 



<3.32) 



and the number o-f groups Cm] is given by 

m = 1 + 1 og 2 N < 3 . 33) 

The groups G 0 and G<log 2 N) represent the power 
content o-f the d.c. component and the -folding -frequency 
component, respectively. The other Groups represent the power 
content o-f the component at n = 1 , n = 2 , n = 3 , etc., plus 
all odd harmonics. 



G 0 = P2<0) 
Gj = P2< 1) 
G2 = P2<2) 
G 3 = P2<4) 
G 4 = P2<8) 
and so on, until 



-»• P2(3) + P2<5) -► ... 

-► P2<6) + P2< 10) -I- ... 

-► P2< 12) -► P2<20) -► . . . 
+ P2<24) P2<40) -*- . . . 



<3.34) 



G<log2N) = P2<N/2) 

The plot o-f Gp versus n gives a discrete spectum 
that is unique -for a certain time sequence C26: p.234], and 
invariant to a time shi-ft o-f the input C23: p. 1063 . 

This spectrum is highly compressed <m points 
instead o-f N /2 -^ 1 points) , and there-fore doesn^'t present all 
characteristics o-f the data as well as the previously de-fined 
spectrum. C23; p.l04] 
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6. The Application o-f the Fast Nal sh Trans-Form 



The application o-f the Walsh Trans-form to 
continuous or discontinuous -functions is easily accomplished 
using the FWT programs provided in Appendix A. The programs 
input a number series (which could represent a sampled 
wave-form) and compute the FWT. With these coe-f-f icients, the 
Power Spectral Density and the Group spectrum can then be 
calculated. The PWT will be applied to the previous examples 
and the results will be compared and analyzed in terms o-f the 
elaborated theory in this chapter. 

a. The FWT o-f a Sinusoid. 

Figure 3.8 gives a normalized graph o-f the 
coe-f-f i c i en ts o-f the FWT o-f one cycle o-f a sinusoid in the 
Walsh interval 0 to 1. The actual magnitudes can be 
determined -from the scale listed on the vertical axis. 

Running either program in Appendix A wi 1 1 give a listing o-f 
the coe-f-f icient magnitudes. 

Note that the maximum coe-f-f icient (normalized 
to 1) is a WAL(n,t) = WAL(l,t). Equation (3.4) o-f the 
chapter reveals there is a SAL( 1 , t) with a sequency o-f 1 at 
this Walsh index. Thus, the one cycle sin wave in the Walsh 
interval 0 to 1 has a distinguishing -feature at a sequency o-f 
1. Perhaps not surprising since it is Known that -frequency is 
a special measure o-f sequency, and the one cycle in this 
interval can be thought o-f as having a frequency of 1. 
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Figure 3.8. F14T o-f 1 Cycle o-f a Sine Function 
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WALSH INDEX (n) 



Also note the other prominent components at 
n = 5, 9, 13, and 29, with sequencies o-f 3, 5, 7, and 15, 

respectively. Our previous calculations with the Walsh series 
(page 71) predicted conspicuous coe-f -f i ci en ts at n = 5 and 



13. 



Applying equations (3.31) to these 
coe-f-f icients results in Figure 3.9, the Power Spectral Density 
(PSD) . The maximum value o-f power is at a sequency of 1, and 
there is no DC power term (sequency of 0) . Harmonic pcwer is 
located at sequencies of 3, 7, and 15. 

These odd number power components combine 
into a single component, Sj, in the Group spectrum. Figure 
3.19. See Equation (3.34) . 

An increase of frequency gives similar 
results. Figures 3.11 to 3.13 show the FWT, PSD, and Group 
Spectrum of three cycles of a sinusoid in the Walsh interval 0 
to 1. The maximum FWT coefficient lies at n = 5 Ca SAL(3,t) 
component with a sequency of 33, with other major ccxnponents 
at n = 1, 9, 13, 21, 25, 29, 53, and 57. All of these 
components are at odd sequencies, so the power will be located 
in the odd number sequencies of the PSD, and the odd numbered 
components will combine into one component of the Group 
Spectrum. Again, see Equations (3.34). 

The FWT, PSD, and Group Spectrum of eight 
cycles of a sine wave in the Walsh interval is analyzed in 
Figures 3.14 to 3.16. The FWT is considerably simpler, with 
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Figure 3.9. PSD o-f 1 Cycle o-f a Sine Function. 
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Figure 3.10. Group Spectrum o-f I Cycle o-f a Sine Function. 
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Figure 3.11. FWT o-f 3 Cycles of a Sine Function. 
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Figure 3.12. PSD o-f 3 Cycles o-f a Sine Function 
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Figure 3.13. Group Spectrum o-f 3 Cycles of a Sine Function. 
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Figure 3.14. FWT of 8 Cycles of a Sine Function. 
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Figure 3.15. PSD o-f 8 Cycles o-f a Sine Function. 
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Figure 3.16. Group Spectrum o-f 8 Cycles o-f a Sine Function. 
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GROUP NUMBER 



only 4 components at PSD sequencies of 8 and 24. The PSD 
shows the peak power at a sequency o-f 8 and a notable harmonic 
at 24. The Group Spectrum -follows -from Equations 3.34. 

A -few conclusions about the trans-form o-f 
a sinusoid could be drawn. The peak coe-f -f i c i en t o-f the FWT 
falls at a Walsh index that translates to a sequency that is 
equal to the sinusoidal frequency. The PSD then has a 
component that is maximum at this sequency. In addition, this 
maximum power component is usually accompanied by notable 
COTiponents at harmonics of this sequency. 

The Group Spectrum, although highly 
compressed and uncomplicated, requires a knowledge of the PSD 
for it to be useful . Note that the Group Spectrum of the 1 
cycle sine wave is equal to the the Group Spectrum of the 3 
cycle wave. An educated guess says that the Group Spectrum 
for all odd numbers of cycles would be the same, that of Gj, 
because only odd numbered ccxnponents <Eq. 3.34) of sequency 
are present. Each component of the Group Spectrum could 
represent a series of sinusoidal frequencies, with no way to 
distinguish the actual frequency of the sinusoid without prior 
knowledge of the PSD. 

b. The FWT of a Rectangular Function. 

Consider the FWT, PSD, and Group Spectrum of a 
4 hertz square wave. Knowing that the Walsh function lends 
itself readily to the reproduction of such a wave, simple 
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results -for graphs o-f the trans-form and related computations 
Mould be expected. 

Indeed this is the case. Figures 3.17 and 
3.18 sh(MM the simplicity o-f represen tat i on -for rectangular 
■functions. The FWT component at n = 7, sequency = 4, is the 
only coe-f-f i ci en t in Figure 3.17. (Note that a 4 hertz square 
wave is WAL(7,t)). All o-f the power lies at a sequency o-f 4 
(Figure 3.18), and in a Group Number 3 (Eqn. 3.34). With a 
proper number o-f samples, any square wave o-f a particular 
■frequency would be represented with similar char ac ter i st i cs . 

A rectangular -function, page 73, composed o-f 
square pulses was expanded in a Walsh Series in this chapter. 
Consider the FWT, PSD, and Group Spectrum o-f the -function. 

Figures 3.19, 3.20 and 3.21 shew the results 
o-f the computation. These graphs are not normalized in 
order to show the similarities with the previously computed 
example on page 73. The FWT shows 2 coe-f -f i c i en ts o-f 0.5 
magnitude, at n = 0 and n = 2, matching the results 
determined be-fore. Since the FWT outputs coe-f -f i c i en ts based 
on positive phasing, this result is correct. 

The power contained in the -function is in the 
DC term and the -first sequency component as shown in Figure 
3.28. These sequency components result in 2 components in the 
Group Spectrum. See Equations 3.34. 
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Figure 3.17. FWT o-f a 4 Cycle Square Wave. 
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Figure 3.18. PSD o-f a 4 Cycle Square Wave. 
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Figure 3.19. 



FWT o-f a 2 Rectangular Pulses. 
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Figure 3.20. PSD of 2 Rectangular Pulses, 
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Figure 3.21. Group Spectrum o-f 2 Rectangular Pulses. 
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The FWT and a Input Time Shi-ft. 

Section 4-c stated that the Discrete Nal sh 



c . 



Trans-form is not invariant to a time shi-ft or phase shi-ft o-f 
the input. Let sin(t) be circularly shi-fted 1t/4 to 
sin<t-1t/4). Now examine Figures 3.22, 3.23, and 3.24 and 
compare with Figures 3.S, 3.9, and 3.19. 

Note the shi-fted FWT coe-f -f i c i en ts pattern is 
more complex. A sizable component now lies at n = 2. This is 
a CAL( 1 , t) CCTnponent with a sequency o-f 1. The unshi-fted FWT 
graph contained only a negligible CAL(l,t> component. 

Double components at n = 5 and 6, 13 and 14, and 29 and 30 
represent sequencies o-f 3, 7, and 15, respectively, the same 
c«nponents that were contained in the unshi-fted sinusoid, but 
the magnitudes o-f the components at the same index Cn] are not 
equal . 

Although the FWT magnitudes are di-f-ferent at 
each index tn] o-f the shi-fted and unshi-fted sinusoids, 
remember that adjacent components combine to -form a PSD 
coe-f-f icient o-f a particular sequency. The components vary in 
a reciprocal manner, so when they are squared and summed, the 
time variant e-f-fect is indeed lessened C23; p.39]. The result 
is a PSD that has sequency components that are very close to 
being time invariant. 

One expectation that was revealed to be true 
was the time invariance o-f the Group Spectrum. Beauchamp C23; 
p.l06] states this spectrum is time invariant to time shifts 



106 



FAST WALSH TRANSFM 
SHIFTED SIN FCN 



I 




Ceut'-T) 3VW 13d 



(D 



<D 

in 



CD 






(T) 



CM 

CD 



m 



Si 

j 



I 



Figure 3.22. FWT o-f a Shi-fted Sinusoid. 
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Figure 3.23. PSD o-f a Shi-fted Sinusoid. 
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Figure 3.24. Group Spectrum o-f a Shi-fted Sinusoid. 
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o-f the input, and examination of Figure 3.23 and its magnitude 
gives an affirmative response to his assertion. 

The rectangular pulse function will positively 
demonstrate the effects of time shift. The FWT , PSD, and 
Group Spectrum of Figure <3. 25a) is shown in Figures 3.19, 
3.20, and 3.21. The FWT and associated spectrums of Figure 
(3.25b) is shown in 3.26, 3.27, and 3.28. 

Quite a difference in the FWT of the shifted 
function can be noticed. The coefficient pattern is much more 
"spread out", and the magnitudes of the prominent coefficients 
at n = 0 and 1 are less in the shifted rectangular function. 

The PSO's are different also. With the power 
being spread among more sequencies in the shifted rectangular 
function, the magnitudes of the pcx*jer coefficients at 0 and 1 
are reduced. 

The Group Spectrums are not alike in this 
case. The two input sequences to the FWT program are made up 
of samples of the original function and the shifted function. 
These sequences are not alike and generate different Group 
Spectrums . 

Has Beauchamp's assertion [23: p. 106] been 
violated? The answer is no. In the sinusoidal shift, it was 
assumed that the function was periodic outside the interval 
0<t<l over which the Walsh functions are defined and over 
which the function was taken. When the sinusoid was shifted, 
the part shifted out of the interval "wrapped around" into the 
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gure 3.25. 



Two Rectangular Pulses, Shi-fted and Unshi-fted. 
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Figure 3 . 26 . FWT o-f Shifted Rectangular Pulses. 
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Figure 3.27. PSD o-f Shifted Rectangular Pulses. 
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Figure 3.28. Group Spectrum o-f Shi-fted Rectangular Pulses. 
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■first part o-f the interval (circular shi-ft), resulting in the 
shi-fted -function having the same power as the unshi-fted 
■function . 

However, the rectangular -function wasn't 
de-fined as a periodic -function, merely two pulses in an 
interval. When the rectangular -function was shi-fted, the part 
that shi-fted outside the interval wasn't wrapped around to the 
■front (not circularly shi-fted) and was lost, decreasing the 
total power -found in the shi-fted -function. The PSD or the 
Group Spectrum could not be the same as the unshi-fted 
■function, since the amount o-f total power is di-f-ferent in each 
o-f the -functions. 

A conclusion is made that the Group Spectrum 
is time invariant to phase shi-fts in the input -for periodic 
circular shi-fted -functions only. 
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IV. THE WALSH TRANSFORM AND DE INTERLEAVING 



The -first two chapters presented some background -for the 
appreciation o-f the role o-f E^ and the basic components o-f a 
typical ESM system. Understanding the tools and purpose o-f the 
trade, though, only gives a small insight into the nature o-f 
probl em. 

There's a lot o-f electromagnetic radiation out there! The 
combat environment presents an overwhelming amount o-f energy 
that has to be collected, analyzed, and sorted. The data 
stream pours into the system -from all angles, with numerous 
-frequencies and amplitudes. Corruptions, disturbances, 
re-flections, and missing pulses make the deinterleaving 
process a most di-f-ficult problem. Today's techniques, 
described in Chapter Two, do a good job, but improvements are 
needed and sought a-fter. 

Unique in many respects, the Walsh Functions, which are 
described in Chapter Three, are examined -for their use-fulness 
in deinterleaving pulse trains into separate chains o-f pulses 
with di-f-ferent pulse repetition intervals. More speci-fical 1 y , 
this chapter attempts to determine i-f there is a PRI/PRF 
recognition -feature in the Walsh Trans-form o-f the data 
presented by an ESM receiver. 
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A. REPRESENTATION OF ESN RECEIVER DATA 



The operating element of the typical deinterleaving 
process is a computer data word (or words) of some format 
appropriate to the system. 

The data word contains values of the initial parameters 
measured by the system, and basically these values are 
examined for matches in previously analyzed pulses. Software 
processing systems make these comparisons and places words of 
similar parameters into cells or histograms, patiently waiting 
until enough pulse information can be found to make an 
identification. Other algorithms process the data word and 
make calculations of PRI's, adding a new sorting parameter 
that is used in the deinterleaving process. 

The investigation of deinterleaving with the FWT will need 
a different symbology for the receiver data, since the FWT 
operates on a series of N numbers, not a data word, to obtain 
the transform. 

Imagine an absolute time line that consists of some point 
in the past, and ends at the present. Consider a received 
pulse as being represented by a vertical line stationed at 
the time of arrival (TOA) of that pulse. The height of the 
line could be made proportional to the amplitude of the 
received pulse. The point on the time line where the pulse 
"stands" is the time of arrival of the pulse. 

Figure 4.1 looks like a messy situation, but there are 
really only four pulse trains of different PRI's on the line. 
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Figure 4.1. TOA representation o-f received pulses. 

PRI's could be calculated by subtracting successive TOA's 
o-f pulses in the same pulse train. 

PRI 1 = T 0 Ai 2 - TOAi 1 
PRI2 = TOA22 - TOA21 
PRI 3 = TOA 32 - TCDA31 
PRI 4 = TOA 42 ~ TOA41 

This procedure is used when other parameters are present 
to identi-fy the pulses. 

Keep in mind that there is a steady stream o-f these pulses 
into the ESM system. In order to make things manageable, the 
number of pulses being examined has to be reduced. 

Lemley 120 : p.l23 calls this reduction mechanism the 
•presorting aperture.' It is of variable width to acccmodate 
high and low density environments and is similar to the time 
slice concept discussed under deinterleaving in Chapter Two. 

It reduces the number of pulses being processed at one time. 

Another way of thinking of this reduction concept is to 
visualize a 'window' on the stream of pulse lines. Only a 
certain number of pulse lines can be present in the window at 
anytime, and deinterleaving operations are carried out on 
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these pulse Ines. The windcxM can be divided into small time 
intervals, which could be analogous to the TOA resolution 
specification o-f a receiver (usually 100—200 nanoseconds). If 
a pulse line is present in the time interval, a number 
proportional to the amplitude is generated for the pulse line 
representation. If not, a zero is used. 

The author prefers to visualize that the receipt of a 
pulse triggers the collection of the next N time intervals for 
analysis and deinterl e'aving. Each interval could contain a 
pulse or not. 

B. CONSTANT AMPLITUDE TOA STRINGS AND THE FWT 

This chapter will work mostly with N = 64 intervals and 
pulse line representations. N must be a power of two to be 
used in the FWT. 

The pulse line representation can be symbolized for 
computer work with the BASIC language DATA statement. Using 
the statement on the 64 pulse lines in the window gives a 
representation that accounts for the two variables that are 
present in the pulse line stream, the amplitude and TOA of the 
pul ses. 

Such a DATA line for Figure 4.1 would be 

DATA 1 , . 3 , . 75 , .3,0, .5, . 75 , 1,0, .3, .5, .3,0, .3, 1, .5,0, .3, . 75 

<4. 1) 

Each number represents the amplitude of a recognized pulse 
that occupies a TOA time interval. 
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Think o-f it in one o-f two ways; 

1. The receiver outputs a measured value o-f a pulse 
amplitude each TOA resolution time interval . I-f no pulse was 
present, it outputs a zero. Or, 

2. The receiver outputs a pulse amplitude only upon 
detection and analysis of a pulse. Time intervals of TOA 
resolution width in between pulse detections can be 
represented with zeroes. 

Either way, the DATA line looks like <4.1). 

Assume for simulation's sake that the period between DATA 
numbers is a known value of time, and the DATA line represents 
the output of an ESM receiver to the preprocessor . The 
numbers indicate amplitude of the received pulses, and their 
distance in intervals between similar pulses is a measure of 
their PRI's. 

Things can be simplified a bit further by considering all 
pulse lines of the TOA string to be of the same amplitude. 

Only one variable would then be present in the string, the 
TOA. 

A DATA line representation of a string of this type is 
DATA 1,1, 1,1, 0,1, 1,1, 0,1, 1,1, 0,1, 1,1, 0,1,1, (4.2) 

In short, this TOA string and its DATA representation is 
merely an indication of the receipt of a pulse or not during 
the resolution time interval. 

A number of Fast Walsh Transforms were run on simulated 
DATA TOA strings, and plots of the coefficients were made. The 
Power Spectral Density and Group Spectrum coefficients were 
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also calculated and plotted. All o-f the Group Spectrum 
plots are not included, as the emphasis o-f this examination 
was placed on the FWT and PSD. Some o-f them are discussed in 
Section C o-f this chapter. 

The plots are grouped by type <FWT or PSD) in Appendix 
D. Each plot is normalized by dividing all coe-f f i c i en ts by 
the maximum coe-f-f icient . Actual values o-f the coe-f -f i c i en ts 
can be determined with the scale on the vertical axis. The 
maximum component's magnitude is shown in parenthesis. 

These plots show the index number and sequency o-f the FWT 
and PSD components, respectively, but only their relative 
magnitude with respect to the maximum component. Keep that in 
mind when examining the plots. Where important, the FWT and 
PSD are replotted with the actual magnitudes on the vertical 
axis. 

These DATA representations simulate single PRI pulse 
trains received by an ESi system. Note the number o-f 
intervals between successive pulses is indicated in the title 
o-f the plot. To calculate a simulated PRI o-f the train, 
simply multiply the number o-f intervals between sucessive 
pulses by the time interval resolution, a known value and a 
-function o-f the receiver. 

The plots are examined by type in the -following sections, 
and conclusions are reached and reported. These single PRI 
representations are the building blocks -for the interleaved 
pulses that are investigated later in the chapter. 



121 



1. Fast Walsh Transforms of the TOA Strings 



An examination of the PWr^s separated them into two 
groups of similar characteristics. 

The first group is composed of FWT's with intervals 
between pulses that are even numbers. Examine each plot 
in Appendix D (page 176 to 197) and note that each even 
interval FWT is symmetr i cal 1 y even about a point midway 
between n = 31 and 32. Also, if the even number is a power of 
two (p.o.t.), the FWT is considerably less ccxnplex. Each 
successive p.o.t. interval is composed of the addition of a 
pair or pair of coefficients to the previous p.o.t. FWT, and 
its magnitude is one half the previous p.o.t. FWT. 

For example, one pair of coefficients is added to TOA 
TAG 2 at n = 31 and 32 (sequency of 16) to form TOA TAG 4. 

Two pairs of coefficients, <n = 15 and 16, n = 47 and 48) are 
added to TOA TAG 4 to make TOA TAG 8, 4 pairs to make TC^ TAG 
16, etc. At N = 64, enough pairs of coefficients have been 
added to have one coefficient for each Cnl. 

The odd numbered intervals do not have symmetrical 
characteristics. In addition, the highest relative magnitude 
Walsh index component or components seem to bear no relation 
to the interval number and thus the PRI of a pulse train. (An 
unfortunate conclusion for PRI recognition.) 

Close scrutiny of both appearance groupings was made 
in an effort to notice a feature that is distinctive to the 
particular similated PRI. There doesn't seem to be any 
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present, either in magnitude or coe-f-f i c i en t distribution, that 
is obvious without a detailed examination. 

0-f initial interest to the author was the magnitude o-f 
the WAL<0,t) ccxnponent. This magnitude seemed to be related 
to the number o-f resolution intervals between the pulses. It 
holds -for most but not all cases. Consider Table 4.1 which 
lists the TAG number and the WAL<0,t) component. 



Table 4.1. Amplitude o-f l*WL<0,t) and TAG No. Inverse 



TOA TAG no. 


NAL<0,t> amplitude 


TAG no. inverse 


1 


1.0 


1.0 


2 


0.5 


0.5 


3 


0 . 34375 


0 . 3333 


4 


0.25 


0.25 


5 


0.203125 


0.20 


6 


0 . 171875 


0. 1447 


7 


0 . 15425 


0 . 1429 


3 


0 . 125 


0 . 125 


9 


0. 125 


0.1111 


10 


0 . 109375 


0 . 1 


11 


0.09375 


0.0909 


12 


0.09375 


0.0833 


13 


0.073125 


0.074923 


14 


0.078125 


0.071429 


15 


0.078125 


0.0447 


16 


0.0425 


0.0425 


19 


0.0425 


0.0524 


22 


0.044875 


0.045455 


25 


0.044875 


0.04 


32 


0.03125 


0.03125 


40 


0.03125 


0.025 


51 


0.03125 


0.019408 


6A 


0.015425 


0.015425 





For 


instance, TOA TAG 


2 has two 


i n terval s 


between 


pu 1 ses , 


and 


the magnitudeo-f WAL<0,t) is 


0.5, or the inverse 


2. TOA 


TAG 


3 has a magnitude 


o-f .34375, 


which is 


close to 
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inverse o-f 3. TOA TAG 4 has a WAL<0,t) amplitude o-f .25, 
which is one divided by -four, the number o-f intervals between 
pulses in the representation. 

I-f there was a relationship between the number o-f 
intervals between the pulses in the DATA line and the 
amplitude, then perhaps this would carry over into the 
interleaved pulse trains composed o-f these single PRI 
representations. Exceptions are evident, though, such as TOA 
TAG 9, which has the same magnitude as TOA TAG 8, although the 
coe-f-f icient distribution is not the same. This implies that 
any recognition -feature will probably depend upon two 
variables, perhaps the amplitude and the sequency o-f the 
coef-f icients. 

2. Power Spectral Densities o-f the TOA Strings 

An examination o-f the PSD's in Appendix D, page 
198 through page 219, will yield similar conclusions. 

PSD's o-f even interval ed TOA strings are symmetr i cal 1 y 
even about a sequency o-f 16. PSD-'s o-f even numbers that are 
also a p.o.t. are considerable less complex, and are one 
-fourth the magnitude o-f the previous p.o.t. PSD (except -for 
TOA TAG 2 to TOA TAG 4) . 

The odd numbered PGDN'S have no symmetry, and sequency 
components that stand out in relative magnitude do not bear a 
relationship to the simulated PRI o-f the string. (Remember 
the number o-f intervals between the TAG'S is a measure o-f the 
PRI o-f the simulated pulse train, when the time interval 
between numbers is Known.) 
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3. General Comments 



The entire purpose o-f this examination has been to 
recognize similarities or relationships between succcessive 
TOA TAG representations in order to recognize the same 
•features <i-f they exist) in an interleaved pulse train DATA 
line representation. 

At this point, the results are very general and 
an optimistic mood regarding any usable -features might belong 
only to the optimist. It seems that the even numbered graphs 
(FWT and PSD ) are the only graphs that show any usable 
features that could be used for recogniton. One shouldn't 
give up all hope, though. 

The author believes it to be possible to develop a 
computer algorithm that would recognize and determine the TOA 
TAG number (and thus a PRI of a simulated pulse train) for a 
given single pulse train DATA line repr esen tat i on . Enough 
features exist in the FWT and PSD that could be used to 
process a particular single PRI represen tat i on with a known 
resolution time interval and output the correct PRI. The 
program would have to determine the PRI by matching 
characteristics (sequency components and magnitudes) of the 
FWT and/or PSD to a data bank composed of characteristics of 
the FWT's and PSD's of single pulse train representations. 

For example, it could first take into account the symmetry of 
the FWT/PSD. If it is symmetrically even, then it isn't an 
odd interval TOA TAG. 
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The algorithm uiould be reasonably complex -for just 
this simple TOA TAG representation o-f a single pulse train o-f 
a particular PRI . But it could be done. 

Since it addresses only single PRI TAG 
representations, the algorithm would only be e-f-fective i-f you 
could convince an opponent to use only one radar at a time 
against the E^1 plat-form. Obviously, this isn't practical. 

More realistically, an examination o-f interleaved 
pulse trains and a determination o-f a PRI recognition -feature, 
i-f one exists, would be in order. This will be done in 
Section D. a-fter discussion o-f the e-f-fects o-f a time shi-ft in 
the TOA TAG string. 

Some limitations o-f the DATA line symbology should be 
stated. Using N = 64 and a resolution between numbers o-f 200 
nanoseconds means the window length is less than 2 
microseconds long. This is not a realistic -figure, since a 
high density o-f 10^ pulses per second gives an average o-f 1 
pulse every microsecond. The use o-f N = 1024 -for an 
interleaved DATA line representation might be -far better in 
the numbers o-f coe-f-f i ci en ts <and more realistic in length o-f 
the window), but unless an identi-f iabl e PRI -feature is 
present, the extra data is exactly that, extra data. 

An additional limitation o-f this representation 
concerns coincident pulses. Pulses arriving at the receiver 
at the same time (-falling in the same resolution time 
interval) are not individually represented in the DATA line 
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symbology. For this reason, one cannot represent a TOA TAG 2 
and any other even numbered TAG. This limitation is a 
problem that will af-fect the FWT and PSD o-f the interleaved 
pulse train representation. 

C. TIME SHIFTED TOA STRINGS 

Chapter Three has stated that the FWT and PSD are not time 
invariant to circular shi-fts o-f the input number series. 
Figures 4.2 through 4.5 are FWT's o-f shi-fted TOA TAGS . 

Compare them with Figures D.4, D.ll, D.13, and D.16. 

The magnitudes o-f the FWT coe-f-f i c i en ts at each Cn] are not 
the same in the shi-fted TAGS. Compare TOA TAG 5 and TOA TAG 5 
SHF 1. However, notice that shi-fted TOA TAG 12 <TOA TAG 12 
SHF 1) has the same magnitude at each WALSH index Cn3 as does 
the unshi-fted TOA TAG, but the components -from n = 32 to n = 

63 are opposite in sign. The even symmetry o-f the even 
numbered FWT about the midpoint has been changed to odd 
symmetry about midpoint. 

This is a very neat and clean distinction, and -further 
emphasizes the distinction between even and odd numbered TOA 
TAGS. It is -further evidence that at least some general 
•features exist which could be made to work in a computer 
algorithm to identi-fy individual PRI"s. 

FWT's o-f odd numbered intervals have no general 
distinctive change between shi-fted plots and standard plots. 

In a -few shi-fted odd numbered plots, there are some particular 
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Figure 4.2. FWT ; TC3A TAG 5 SHIFT 1 
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WALSH INDEX (n) 
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magnitudes o-f index Cnl that remain the same or increase 
proportionate! y , but others do not. I -f any real relationship 
between the shifted and unshi-fted odd numbered FWT plots 
exist, then it will have to be revealed with a much more 
detailed study. 

Quite simply, odd numbered shifted FWT plots show no major 
distinctive features that are readily apparent. Even numbered 
shifted FWT plots follow the sign change rule for components 
at n = 32 to 63, but the magnitudes of the components 
of both plots are the same. 

Although the odd numbered PSD plots change in the 
amplitude of the sequency components from unshifted to 
shifted, the even numbered PSD's are the same in both shifted 
and unshifted plots. Indeed, they are the same in both 
magnitude and location of sequency components. Compare 
Figures 4.6 to 4.9 on the following pages with Figures D.26, 
D.33, D.35, and D.38 in Appendix D. 

Thus, the even numbered PSD^’s of this DATA line 
representation join the time invariant Group Spectrums 
in being immune to the effects of a circular time shift. 
Compare the shifted Group Spectrums Figures 4.10 to 4.13, with 
the unshifted Group Spectrums in Appendix E, Figures E. 1 
through E.4. 

Again the point is made that these general features could 
be recognized in a PRI recognition algorithm. At the very 
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least, it would eliminate the odd or even choices o-f time 
intervals between numbers o-f the single PRI representation. 



D. INTERLEAVED TOA STRINGS AND THE FWT 

The receipt of radar pulses from several different 
emitters results in individual pulse trains that are 
interleaved into a data stream from the receiver. Separation 
or deinterleaving of these pulse trains from the data stream 
must be done to identify the individual emitters. This 
section examines the plots of the FWT and PSD coefficients of 
a DATA line represen tation of interleaved pulse trains 
composed of the TOA TAG repr esen tat i on in Section B. of this 
chapter . 

It is desired to answer the following questions; 

1. Are the components of the FWT and PSD of an 
interleaved represen tat i on identifiable as belonging to a 
particular PRI? 

2. Is there a relationship between the magnitude or 
sequency of the components of an interleaved plot with that 
of plots of single PRI trains? 

3. Are interleaved plots made by simply adding single PRI 
pi ots? 

Begin by considering interleaved TOA TAGS of Table 4.2. 
The DATA line representati on for TOA TAG 5 wi th N = 64 is 



DATA 1,0, 0,0, 0,1, 0,0, 0,0, 1,0, 0,0, 0,1, 0,8, 0,0,1 
0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0,0 
0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0,0 



0 , 0 , 0 , 0 , 1 ,0 
0 ,0 , 1 , 0 , 0 ,0 



<4. 3) 



and for TOA TAG 8 




0 , 0 , 0 , 1 , 0 ,0 
1 , 0 , 0 , 0 , 0 ,0 



<4. 4) 
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Now, interleave the two, 

DATA 1,0, 0,0, 0,1, 0,0, 1,0, 1,0, 0,0, 0,1, 1,0, 0,0, 1,0, 0,0, 1,1,0, 
0 , 0 , 0 , 1 , 0 , 1 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 1 , 0 , 1 , 0 , 0 , 0 , 
0,1, 1,0, 0,0, 1,0, 0,0 <4. 5) 

Notice that the first pulse of each TAG is coincident and 
both are represented by the first pulse of the interleaved 
DATA line. 

The FWT of TAGS 5 and 3 are shown in Figures 4.14 and 
D.7. TOA TAG 5 components are emphasized by setting the 
WAL<0,t) components equal to 0 before plotting. The plotting 
routine then normalizes to a smaller magnitude component and 
increases the relative magnitude of the other components. 

Examine the interleaved FWT, Figure 4.15. Note sizable 
ccmponents at 9, 22, 25, 31, 32, 37, 45, 48, 50, 60, and 63. 
Now examine TOA TAG 5 with WAL<0,t) equal to 0. Note 
emphasized components at 25, 37, 44, 45, 50, 51, and 52. 

There''s some commonality between the plots. Components at 
25, 37, 45, 50, and 52 are obviously associated with TOA TAG 
5. So where is TO^ TAG 3? A closer look reveals its 
influence on the interleaved FWT. Note identical relative 
magnitude pairs at n = 15 and 16 and n = 31 and 32. These are 
generated by the coefficient pairs of these Cn] from TOA TAG 
8. Note also that n = 48 of the interleaved FWT is higher in 
relative magnitude because there is a small n = 48 component 
of TOA TAG 5. 

How would one use this interleaved FWT to deinter leave 
pulse trains? A direct computation method using a correlation 
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algorithm is considered in the next section. The author 
suggests a "matching" algorithm could be written that would 
scan the interleaved FWT, and note relatively high magnitudes 
at certain Cn] above some threshold. These Cn] could be 
matched to a data bank o-f important Cn] components o-f single 
PRI FWT plots, along with any special char ac ter i st i cs o-f the 
interleaved FWT (identical components at successive Cn] , -for 
example). Promising matches between the interleaved 
parameters and single PRI parameters would result in a 
decision or choice o-f which single PRI trains make up the 
interleaved pulse train. 

Scrutiny o-f the PSD plots generates similar conclusions. 
Figure 4.16 shows the interleaved PSD plot. Components o-f 
relative magnitude o-f 0.3 or greater are 3, 13, 16, 19, 22, 

24, 25, 26, and 32. 

Now note the PSD's o-f T(^ TAG 5 and 8 (Figures 4.17 and 
D.29. TOA TAG 5 contributes components at sequencies o-f 13, 
22, 25, and 26. TOA TAG 8 contributes components at 8, 16, 24, 
and 32. 

Be-fore moving on to another interleaved case, a discussion 
o-f the magnitudes o-f the interleaved components (both FWT and 
PSD) is in order. A de-finite relationship is not apparent 
between components o-f individual plots and interleaved plots. 

The author can say that one's intuition about components 
o-f interleaved plots generally holds true. For instance, a 
large positive FWT component of one TAG and a negative FWT 
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component o-f another TAG at the same Cn] will usually generate 
a smaller canponent in a FWT component o-f the interleaved 
TAG'S. The result is not a simple addition, but the actual 
magnitude o-f the interleaved component is usually smaller than 
the added values o-f the single ccMnponents. No -fast and easy 
rul e appi ies. 

In a matching algorithm, the magnitude would have only a 
comparative bene-fit. The actual promise of the algorithm lies 
in matching component locations (either Cn3 or sequency) to a 
data bank of single PRI ccMnponent locations, with magnitude 
and other special character i sties as a decision making aid. 

One more interleaved case is examined. This one is 
different because there are no coincident pulses. The shift 
in TOA TAG 5 causes it to interleave perfectly with TOA TAG 
19. 

TOA TAG 5 SHF 1 (Figure 4.2 and 4.6) DATA line; 

DATA 0,1, 8, 0,0, 0,1, 0,0, 0,0, 1,0, 0,0, 0,1, 0,9, 0,0, 1,0, 9, 0,9, 1,0, 
0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 
0,1, 0,0, 0,0, 1,0,0 (4.5) 

TOA TAG 19 (Figure 4.18 and 4.19) DATA line; 

DAT A 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 1,0, 0,0, 0,0, 0,0,0, 
0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 
0,0, 1,0, 0,0, 0,0,0 (4.6) 

And the interleaved DATA line: 

DATA 1,1, 0,0, 0,0, 1,0, 0,0, 0,1, 0,0, 0,0, 1,0, 0,1, 8, 1,0, 0,0, 0,1,0, 
0,0, 0,1, 0,0, 0,0, 1,0, 1,0, 0,1, 0,0, 0,0, 1,0, 0,0, 0,1, 0,0,0, 
0,1, 1,0, 0,0, 1,0,0 (4.5) 

Figure 4.20 and 4.21 show the FWT and PSD of the 

interleaved TAG'S. Interleaved components of both single PRI 
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Figure. 4.28. FWT ; TOA TAG 5 AND 19 INTERLEAVED, 
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TAGS are present, and the same process as be-fore is used to 
identi-fy them. A new discovery o-f this examination reveals 
that some -feature PSD components o-f TOA TAG 5 SHF 1 are passed 
directly to the interleaved PSD. Note the magnitude of 
sequencies 11, 16, and 21 of the interleaved PSD is the same 

as the PSD of TOA TAG 5 SHF 1, and that the magnitudes of TOA 
TAG 19 PSD at these sequencies is zero. This is the first 
evidence of two single train represen tati ons generating an 
interleaved components magnitude by merely adding the simple 
PRI components. A more closer study could be done to see if 
it is merely coincidence. The real question is," Did the fact 
that there are no coincident pulses in the interleaved 
representation allow me to add certain single PRI components 
to form an interleaved component?" 

A partial answer to the question probably is evident since 
only certain components of the two single PRI PSD plots have 
this feature of simple addition. 

E. CORRELATION THE PSD COEFFICIENTS 

Throughout this chapter the emphasis has been placed on 
the location and relative magnitudes of components located in 
the Malsh index domain or the sequency domain. Most of the 
plots show the magnitude of the components relative to the 
maximum component magnitude of that particular plot. It was 
stated several times that the actual magnitudes of the 
individual PRI (TOA TAG) components didn^'t seem to have a 
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recognizable relationship with the magnitudes o-f the 
components o-f the interleaved plots. 

It was thought that the magnitudes might be related, and a 
correlation o-f the interleaved ccMnponent magnitudes with those 
o-f the individual PRI magnitudes might show their presence. 

The PSD coe-f-f i c i en ts were chosen because they are de-fined -from 
the FWT coe-f-f i c ien ts , and exhibit the same characteristics and 
in-formation about the trans-form in terms o-f symmetry and 
relative magnitude components. 

Without any mathematical -fan-fare, the correlation 
coe-f-f icient P is presented. 

P = E < <X-X/)/dx • <Y-Y/)/<?y > <4. 

where 

X,Y are random variables 

X/, Y^ are the mean of the variables 

variances o-f the variables 

The correlation coe-f-f icient expresses the degree to which 
two random variables are correlated without regard to the 
magnitude o-f either one. C30: p.83I It is a normalized 
quantity and will be a value between -1 and -► 1 . 

In this application, X are the PSD coe-f -f i c i en ts o-f the 
interleaved situation and Y are the PSD coe-f-f ic ien ts o-f the 
individual TAG situation. 

A correlation BASIC language routine was written 
speci-f ical 1 y -for the 32 PSD components o-f the interleaved 
and individual TOA TAGS. This program was adapted -frtxn a 
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correlation routine in the IMSL library contained in the Naval 
Postgraduate School 3933 computer. It was not adapted -for 

general application but spec i -f i cal 1 y -for this situation. For 
this reason it is not provided in Appendix A. 

The correlation routine merely reads in the component 
values o-f the individual PRI TAG representations and the 
interleaved representation generated by the FWT/PSD/Group 
Spectrum program provided in Appendix A. It correlated the 
two, based upon Equation (4.5), and outputed a single number, 
the correlation coe-f-f icient . 

The PSD coe-f-f icients o-f the interleaved situation o-f TOA 
TAG 5 and TOA TAG 1? were used (Figure E.5 in Appendix E) . 

This interleaved situation is di-f-ferent -from that o-f page 148 
because TOA TAG 5 wasn't shi-fted be-fore interleaving, causing 
the -first pulse o-f TOA TAG 5 and 19 to be coincident. With 
these coe-f -f i c i en ts , the values o-f the individual TOA TAG 
components -from 3 to 64 were then correl ated with the 
interleaved values using the correlation routine. 

Table 4.2 gives the TOA TAG number, and the correlation 
coe-f-f icient generated. 

Examination o-f the correlation coe-f-f icients shows strong 
correlation with TOA TAG 5. Great! It is desired to have a 
high correlation value with TOA TAG'S 5 and 19, since the 
interleaved representation was composed o-f these two TAG 
represen tat i ons . The TOA TAG 5 correlation coe-f-f icients -fits 
the need nicely. 
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Table 4.2. Correlation o-f Interleaved PSD coe-f -f i c i en ts with 

individual TOA TAGS 



TOA TAG number 



Correlation Coefficient 



3 


0.665 


4 


0.234 


5 


0.963 


6 


0.582 


7 


0.305 


8 


0.030 


9 


0 . 166 


10 


0.709 


1 1 


0.299 


12 


0 . 183 


13 


0.274 


14 


0. 191 


15 


0.365 


16 


-0.065 


19 


0. 140 


22 


0.330 


25 


-0.028 


32 


-0.055 


40 


0. 106 


51 


0.059 


64 


-0.722 



On the down side, the interleaved PSD coefficients do not 
show a high value of correlation with TOA TAG 19. The initial 
success of TOA TAG 5 certainly indicates that this is an area 
for further investigation, however. 

The author suggests that the correlation procedure could 
be used in an effective deinterleaving algorithm if it was 
applied to the sequential method of extracting PRI's. Once the 
initial correlation routine is carried out and a high value of 
correlation coefficient is identified and associated with an 
individual TOA TAG string, then the algorithm could return to 
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the interleaved DATA stream representation and subtract (set 
equal to zero) those TOA TAG'S associated with the determined 
PRI . 

Now repeat the entire process. Take the FWT and PSD o-f 
the remaining TOA TAG r epresen tat i on , correlate the PSD 
coe-f-f icients with the data bank o-f individual TOA TAG PSD 
coe+-f ici ents, and see if a new TOA TAG has been identi-fied 
with another high value o-f correlation coefficient. If it 
has, then strip these TOA TAG'S from the original string. 

The process could be repeated until calculation of the PSD 
yields no further useful coefficients, or until no correlation 
IS found. The stripped TAG's represent the individual pulse 
trains that make up the interleaved pulse string. 

This suggested process is similar to the way a 
preprocessor sorts the incoming data stream. After 
identifying particular pulses in the stream, the preprocessor 
strips these pulses from the stream, simplifying the follow-on 
process, which continues to identify < dein ter 1 eave) individual 
pulse trains of the remaining stream. Each identified train 
is then removed from the incoming data from the receiver. 

A problem that plagues current deinterleaving techniques, 
that of coincident or missing pulses, would probably plague 
this correlation deinterleaving procedure. 

The TOA TAG representation does not indicate coincident 
pulses separately. Coincident pulses are represented by a 
single TOA TAG. Only the receiver knows which one it detected 
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and represented as a received pulse. There-fore, the stripping 
o-f this TAG generated when two pulses had arrived at the 
receiver at the same time inadvertently strips a TOA TAG -from 
the remaining pulse representation . 

I-f the PSD coefficients of the remaining string with the 
inadvertently stripped TAG were now calculated, one would find 
they are different than those calculated from the remaining 
string that included the inadvertently stripped TAG. 

These PSD coefficients, those calculated frcm the missing 
TAG string, cannot correlate higher with the individual TOA 
TAG strings (used by the correlation routine) than the 
coefficients of the complete remaining string. 

The individual TAG PSD coefficients were calculated from 
T(^ TAG strings with no missing pulses. So, coefficients from 
strings with missing pulses will not correlate as well as 
coefficients calculated from complete strings. 

A possible solution could be the use of a large number of 
TOA TAG'S in the interleaved represen tat i on , say 1024 or 2043. 
In other words, expand the "window" or aperture on the data 
stream. This would allow the interleaved representation to be 
more realistic and comparable to the receiver's own TOA 
resolution. More coefficients might help the correlation 
procedure and lessen the effect of missing and coincident 
pu 1 ses . 
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Kf. CONCLUSIONS AND RECOMMENDATIONS 



Be-fore discussing the conclusions o-f the research e-f-fort 
and making recommendations -for -follow-on work, a review o-f the 
investigation steps that were taken by the author and 
translated into the substance o-f this thesis will be 
bene-ficial . 

For background, the importance o-f the ESM e-f-fort and its 
place in Electronic War-fare was studied. A typical ESM 
system was examined, and the problem o-f deinterleaving the 
pulse stream output -from an ESM receiver was considered, along 
with current techniques -for deinterleaving already in use. 

Having no previous idea o-f even the existence o-f the Walsh 
-functions set the stage -for the study o-f these unique and 
interesting -functions. As the study proceeded through the 
de-finition o-f the Walsh Trans-form, the constant comparisons o-f 
the Walsh Trans-form properties with the -familiar properties o-f 
the Fourier Trans-form brought review o-f this area. 

Re-ference 23 provided the FWT program that became the tool 
o-f the investigation. Adaptation o-f the program to the BASIC 
language and a gathering o-f evidence that supported belie-f in 
its output took a considerable time. Searches o-f the 
literature -for examples o-f FWT's o-f -familiar waveshapes 
-finally provided the proo-f o-f the program. 

With the guidance o-f the thesis advisor, the 
representation o-f the receiver data stream with time o-f 
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arrival tags was accepted, and the FWT program was applied to 
many individual tag strings. The purpose was to determine if 
a TOA TAG string representing a particular PRI pulse train 
could be recognized by certain features of the Walsh Transform 
or the Power Spectral Density. The FWT, PSD, and Group 
Spectrum of interleaved TOA TAGS were then closely examined 
for PRI recognition features. Finally, a correlation 
algorithm was used to determine the degree of correlation 
between interleaved PSD coeffcients and single pulse train PSD 
coefficients. An examination of the correlation coefficients 
was done to see if a PRI could be determined by noting the 
magnitude of the coefficient. 

A. CONCLUSIONS 

One result of the effort was the demonstration of the 
ease of computation of the Walsh Transform and Power Spectral 
Density of a number series. The algorithm is fast and 
efficient, and easily adapted to a BASIC language 
implementation. This allows the computation of the transform 
and calculation of the PSD coeffcients on a microcomputer , 
easily used and readily available in the work space. 

When applied to a TOA TAG representation of a single PRI 
pulse train, the transform shows features that can be used to 
recognize this pulse train when compared to the transforms of 
other PRI pulse trains. Symmetry of the coefficient 
distribution, and the relative magnitude of certain Walsh 
index or sequency components are the features that allc^Ai this 
identification . 
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It was shown in Chapter 4 that each single PRI pulse 
train's Walsh transform has unique and distinctive components 
that distinguish it -from other transforms of single PRI 
representations. No single feature, however, such as a 
relative maximum sequency component, can be connected with the 
PRI of the pulse train. Each particular PRI transform simply 
has individual features that identify it, and only certain 
characteristics, such as symmetry, are cofranon between groups 
of individual pulse train transforms. 

The Walsh transform and PSD's of interleaved TOA TAG 
representations have features that can be identified as 
belonging to the transforms or PSD's of component individual 
pulse trains. However, these features are not readily 
apparent and one could not identify or associate them with a 
particular single PRI transform without prior Knowledge of 
single PRI transform features. Like the transforms of single 
PRI pulse trains, the interleaved transform or PSD does not 
have any features that can identify component PRI's from the 
transform or PSD alone. 

This lack of component PRI features in the transform or 
PSD effectively prevents its use in a practical deinterleaving 
algorithm. At the outset of the research effort, it was hoped 
that the transform of a representation of the pulse stream 
from a E^ receiver would identify the individual PRI's of the 
component pulse trains. Neither the transform or the PSD does 
this, although there are cannon features between the 
interleaved transform/PSD and component tr ansf orm/PSD . 
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The application o-f a correlation routine between the 
the interleaved PSD coe-f-f icients and the individual TOA 
TAG PSD coefficients was seen to be effective in identifying 
one of the component PRI's. It did not identify all 
component PRI's, and any algorithm that uses this 
approach will have to go through additional steps to do so. 
These required additional steps quickly negate the speed of 
computation advantage that the Walsh transform exhibits, and 
the approach becomes quite similar to current sequential PRI 
extraction techniques. 

In summary, the transform of the interleaved TOA TAG 
representation failed to produce a simple PRI recognition 
feature. Recognition features are present, however, in the 
interleaved transform and PSD that would allow the 
identification of component pulse trains in an interleaved 
representation if a data bank of individual PRI transf orm/PSD 
features was available. 

B . RECOMMEWDAT I ONS 

It is recommended that an effort at writing a "matching 
feature" deinterleaving program be made. The program would 
use a data bank of parameters that describe the interesting 
properties and features of the Walsh transform and/or PSD of a 
TOA TAG represen tat i on of a single PRI pulse train. 

The data bank would be the measure of attaining a 
successful deinterleaving program. Properties such as even 
symmetry, double components of power of two interval 
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representations, and distinctive index or sequency components 
o-f the odd numbered TOA TAG representations would be 
included. 

The algorithm would scan the FWT/PSD o-f the interleaved 
TOA TAG representation -for recognizable -features, then match 
these -features to the data bank. Appropriate tolerances could 
be included to be used in decisions about whether a single PRI 
is a component o-f the interleaved train. 

Additional study could possibly yield the optimal -features 
to include in the bank. For certain, a closer look at single 
PRI representations must be done to determine which components 
are shared among them. The author -firmly believes that there 
will always be a FWT/PSD -feature o-f a single PRI 
represen tati on , whether a particular index or sequency 
component, or a symmetrical charac ter i st i c , that can be used 
to distinguish it -from other single PRI representations. 

Closer study o-f whether there is a relationship between 
interleaved and individual index or sequency magnitudes should 
be COTipleted. It is a -feature that could possibly be used in 
a matching algorithm, and it certainly can be used in direct 
compuation methods, such as a least mean squares or 
correlation routine that might generate an individual PRI 
indication . 

The quite tentative but interesting result -found in the 
last section o-f Chapter 4 should be explored completely. All 
coe-f-f icients should be involved, including the Group Spectrum 
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€ 



m 




coe-f -f i c i en ts 



. In this thesis, the Group Spectrum c 
have mainly been used to show the time invariance n 
their character, but they might o-f-fer some interest 
computation -features. 



oe-f -f i c i en ts 
ature o-f 
ing direct 
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APPENDIX A 



FAST NALSH TRANSFORM PROGRAMS 

The follcftaing pages contains listings o-f the 
Walsh Transform programs used in this thesis. The first 
program is written for an IBM Personal Computer, and the 
second was written for the HP-85. 

The FWT subroutine was adapted from a FORTRAN program in 
■Walsh Functions and their Applications," by K.G. Beauchamp. 

The PSD and Group Spectrum subroutines were written by the 
author . 

Both programs were run on the respective computer BASIC 
interpreter (uncompiled). A compiled version would probably 
run faster, although a 128 number FWT ran on the IBM PC in 
about 4 seconds using only the interpreter. 
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10 '' This program computes the Fast Walsh Tran-form , the 
Power Spectral Density, and the Group Spectrum 

20 ' Coe-f-f icients o-f a number series. 

30 ' The FWT subroutine is a BASIC adaptation o-f a FORTRAN 
program listed 40 ' i n “ Walsh Functions and Their 
Applications," by K.G. Beauchamp. 

45 ' The PSD and Group Spectrum subroutines were written by 
the author. 



50 

60 

65 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

230 

235 

240 

250 

260 

270 

280 

290 



' This program runs on a IBM Personal Computer. The 
printer used was a NEC 8023. 

DIM X< 100) ,Y< 64) , P< 64) , G< 1 1) ,F< 64) 'dimension arrays 
LPRINT CHRS<27) ;"L" ; "008" 'set le-ft margin 

N = 64 ' number o-f input numbers 

FOR I = 1 TO N ' read in data 

'READ X<I) '-from the DATA line, or 

XU) = SIN< 13fI/N5f2X3. 1416) 'calculate your own 
NEXT ' I 

G0SU6 310 'calculate FWT o-f input series 

FOR 8 = 1 TO N 'divide all coe-f-fs by N and downshft 
F<B-1) = X(B)/N ' this array contains the FWT coe-f-fs 
NEXT 'B 

'calculate power spectrum coe-f-f i ci en ts 
' calculate group spectrum 
'output the FWT coe-f-fs 
'output the PSD coe-f-fs 
'output the group coe-f-fs 



GOSUB 

GOSUB 

GOSUB 

GOSUB 

GOSUB 

END 



600 

720 

840 

950 

1047 



SUBROUTINE FHT<N,X,Y) 



This routine per-forms a FWT o-f an input series in array 
X. The array Y is used -for working space. 

The dimensions o-f X and Y must be a power o-f 2. 

The results o-f the FWT are in sequency order , positive 
phasing, and in array X. The subroutine uses a Hadamard 



trans-form . 

300 ' 

310 N2 = N/2 

320 M = 6 ' M equals the base 2 logarithm o-f N 

330 FOR L = 1 TO M 



340 


NY = 0 








350 


N2 = 2 


'^<L 


-1) 




360 


N2I = 


2 % 


NZ 




370 


N2N = 1 


N/'NZI 




380 


FOR I 


= 1 


TO 


NZN 


390 


NX = 


NY 


+ 


1 


400 


NY = 


NY 


* 


NZ 


410 


JS = 


U 


-1) 


X NZI 


420 


JD = 


JS 


-I- 


NZI + 


430 


FOR 


J = 


NX 


TO NY 


440 


JS 




JS 


+ 1 
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450 J2 = J + N2 

460 Y<JS> = X<J) + X<J2> 

470 JD = JD - 1 

480 Y(JD> = X<J)- X<J2> 

490 NEXT ' J 

500 NEXT ' I 
510 FOR B = 1 TO N 
520 X(B> = Y<B) 

530 NEXT ' B 
540 NEXT 'L 
550 RETURN 

555 ' SUBROUTINE PSD<N) 

560 ' 

570 ' 

580 ' This routine calculates the Walsh Power Spectral 
Density Coe-f-f icients. 

590 ' 

600 P<0> = <F<0))''2 'first PSD coefficient 

610 K=1 

620 FOR I = 1 TO <N/2)-l 'coefficients 1 to N/2-1 
630 P<I)= <F<K))''2 + <F<K+1))''2 

640 K = K + 2 

650 NEXT 'I 

660 P<N/2>=<F<N-1)>''2 'last PSD coefficient 

670 RETURN 
680 ' 

690 ' SUBROUTINE GROUP SPECTRUM<M) 

700 ' 

710 ' This subroutine calculates the group spectrum 
coefficients. 

720 D=1 

730 G(0) = P<0)^2 
740 FOR B = 0 TO M-1 

750 FOR C =<2'"B) TO N/2-1 STEP (2X2'‘B) 

760 G<D) = G<D) + P<C)^2 
770 NEXT 'C 
780 D = D+1 
790 NEXT 'B 
800 G<M> = P<N/2)''2 
810 RETURN 
820 ' 

830 ' This subroutine outputs the Walsh Transform 
Coef f icients. 

840 ' 

841 FOR C = 1 TO 5: LPRINT; NEXT ' 1 inch top margin 

842 LPRINT" FWT COEFFICIENTS FOR 1 CYCLE SINE 

FUNCTION" 

843 LPRINT; LPRINT 'skip two lines 

850 LPRINT "Walsh Function" TAB<20) "Coeff." TAB< 35) "Wal sh 
Function" TAB(55) "Coeff." 
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I 




LPRINT; NEXT 
'■form -feed 



'skip 16 lines 

xx" 'print pg. no 



860 LPRINT STRING^C 15,45) TABC20) STRINGS< 6 ,45) TAB<35) 
STRING*< 15,45) TAB<55) STRING$< 6 , 45) 

870 FOR B = 1 TO 32 'output the coe-f-fs 

880 LPRINT “WAL<" ;B-1;" ,t) “ TAB< 16) ; LPRINT USING 

;F<B-1) ; : LPRINT" " TAB<35) 

890 LPRINT "WAL<" ;B+31;“ ,t) ■ TAB<52) : LPRINT USING 

;F<B■^31) 

900 NEXT 'B 

910 FOR C = 1 TO 16 

911 LPRINT" 

912 LPRINT CHR*<12) 

920 RETURN 

930 ' 

940 ' Output the PSD coe-f -f i c i en ts . 

950 ' 

951 FOR C = 1 TO 5; LPRINT: NEXT 

970 LPRINT " PSD COEFFICIENTS FOR A 1 CYCLE SINE 

FUNCTION" 

980 LPRINT : LPRINT 
990 FOR B = 0 TO N/2 

1000 LPRINT "P<"; B ; ")=" LPRINT USING P< B) 

1020 NEXT 'B 

1025 LPRINT: LPRINT 

1040 RETURN 

1044 ' 

1045 ' Output the group coef -f i c i en ts 

1046 ' 

1047 LPRINT- GROUP SPECTRUM COEFFICIENTS OF A 1 CYCLE SINE 

FUNCTION" 

1048 LPRINT: LPRINT 
1050 FOR B = 0 TO M 

1060 LPRINT -G<-;B;-)= ";: LPRINT USING " ; G< B) 

1070 NEXT 'B 

1071 FOR C = 1 TO 7: LPRINT: NEXT 

1072 LPRINT - XX- 

1073 LPRINT CHR*< 12) '-form -feed 
1080 RETURN 
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10 ! This program ccmputes the Fast Walsh Tran-form , the 
Power Spectral Density, and the Group Spectrum 

20 ! Coefficients of a number series. 

30 ! The FWT subroutine is a BASIC adaptation of a FORTRAN 
program listed 40 ! in " Walsh Functions and Their 
Applications,” by K.G. Beauchamp. 

45 ! The PSD and Group Spectrum subroutines were written by 
the author. 



50 ! This program 
coef f i c i en ts 
DIM X< 100> ,Y<64> 
PRINTER IS 2 
N = 44 ! 

FOR I = 1 TO N 
!READ X<I> 

XU) 

NEXT 



40 

45 

70 

80 

90 

100 

110 

120 

130 

140 

150 

155 

154 

140 

145 

170 

175 

180 

190 

200 

205 

204 

207 

208 

209 

210 
211 
212 

213 

214 

215 
214 

217 

218 

219 

220 
221 
222 
225 

230 

231 



runs on HP 85 Computer, and writes the 
to data files named “ FWT" , " PSD" , and “GR 
, P< 44) , G< 1 1) , F< 44) !dimension arrays 

! select printer for output of print 
number of input numbers 

! read in data 

! f rom the DATA line, or 

= SIN< na/N*23f3. 1414) Icalculate your own 
I 



Icalculate FWT of input series 
! divide all coeffs by N 
! this array contains the FWT coeffs 



GOSUB 310 
FOR B = 1 TO N 
X(B) = X<B)/N 
NEXT B 
BEEP @ CLEAR 
DISP ”FWT CAL OVER" 

GOSUB 400 Icalculate 

DISP "PSD CAL OVER" 

GOSUB 720 ! calculate 

DISP "GRP COEFF CAL OVER" 

GOSUB 840 ! output the FWT coeffs 

GOSUB 950 ! output the PSD coeffs 

GOSUB 1047 ! output the group coeffs 



power spectrum coefficients 
group spectrum 



CREATE 
ASSIGN 
FOR B = 1 
PRINT# 1; 
NEXT B 
ASSIGN# 1 
CREATE 
ASSIGN #1 
FOR B = 0 
PRINT# 1; 
NEXT B 
ASSIGN# 1 
CREATE 
ASSIGN 
FOR B = 0 
PRINT# 1; 
NEXT B 
ASSIGN# 1 



"FWT" ,30 
#1 TO "FWT' 
TO 44 
X<B) 

TO * 

PSD" , 20 
TO "PSD' 
TO N/2 
P<B) 



I creates file 
opens file 



closes file 



TO X 
"GRP" , 
#1 TO 
TO M 
G<B) 

TO * 



20 

GRP' 



BEEP 

END 

I 



e DISP "PROGRf^ OVER' 
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SUBROUTINE FHT<N,X,Y> 



235 ! 

240 ! 

250 ! 

260 ! This routine per-forms a FWT o-f an input series in array 
X. The array Y is used -for working space. 

270 ! The dimensions of X and Y must be a power of 2. 

280 ! The results of the FWT are in sequency order, positive 
290 ! phasing, and in array X. The subroutine uses a Hadamard 
transform. 

300 ! 

310 N2 = N/2 

320 M = 6 ! M equals the base 2 logarithm of N 

330 FOR L = 1 TO M 

340 N1 = 0 

350 N3 = 2^<L-1> 

360 N4 = 2 X N3 

370 N5 = N/N4 

380 FOR I = 1 TO N5 

390 N6 = N1 + 1 

400 N1 = N1 + N3 

410 J1 = <I-1) * N4 

420 J3 = J1 + N4 + 1 

430 FOR J = N6 TO N1 

440 J1 = J1 + 1 

450 J2 = J + N2 

460 Y<J1) = X<J> + X<J2) 

470 J3 = J3 - 1 

480 Y<J3) = X<J>- X<J2) 

490 NEXT J 

500 NEXT I 

510 FOR B = 1 TO N 

520 X<B> = Y<B> 

530 NEXT B 
540 NEXT L 
550 RETURN 

555 ! SUBROUTINE PSD<N) 



580 ! This routine calculates the Walsh Power Spectral 
Density Coefficients. 

590 ! 



560 ! 
570 ! 



600 P<0) = <F<0>>^2 
610 K=1 

620 FOR I = 1 TO <N/2)-l 
630 P<I)= <F<K))^2 + 

640 K = K + 2 

650 NEXT ! I 
660 P<N/2)=<F<N-1) > -^2 

670 RETURN 

671 ! 

672 • 



! coefficients 1 to N/2-1 
<F(K+1>> ^2 



•.last PSD coefficient 



[first PSD coefficient 
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SUBROUTINE GROUP SPECTRUM<M) 



680 
690 
700 

710 ! This subroutine calculates the group spectrum 
coefficients. 

720 D=1 

730 G<0) = P<0)'‘2 
740 FOR B = 0 TO M-1 

750 FOR C =<2"B) TO N/2- 1 STEP <2?52'‘B) 

760 G<D) = G<D) + P<C)'‘2 

770 NEXT !C 

780 D = D+1 

790 NEXT B 

800 G<M) = P<N/2)^2 

810 RETURN 

820 I 

830 ! This subroutine outputs the Nalsh Transform 
Coefficients. 

840 I 

850 PRINT "Walsh Function" ;TAB<20) ; "Coefficient" @ PRINT 
870 FOR B = 1 TO 64 ! output the coeffs 

880 PRINT "WAL<" ;B-1;“ , t) " ;TAB<20) ; F(B) 

900 NEXT B 
920 RETURN 
930 ! 

940 ! Output the PSD coefficients. 

950 ! 

960 PRINT 

970 PRINT "POWER SPECTRUM COEFFS" @ PRINT 

990 FOR B = 0 TO N/2 

1000 PRINT "P<"; B ;")=* ; P<B) 

1020 NEXT B 
1040 RETURN 

1044 ! 

1045 ! Output the group coefficients 

1046 ! 

1048 PRINTS PRINT 

1049 PRINT "GROUP SPECTRUM COEFFICIENTS" S PRINT 

1050 FOR B = 0 TO M 

1060 PRINT "G<";B; •)= SPRINT USING G< B) 

1070 NEXT B 
1080 RETURN 

1085 ! This is a typical DATA line of numbers representing 
TOA TAGS 

1090 '.DATA 1,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0,1, 
0 ,0 ,0 ,0 ,0 ,0 ,0 , 1 "etc. for N = 64 
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APPENDIX B 



ORTHOGONALITY 

The concept o-f orthogonality in a set o-f -functions is 
important because only orthogonal sets o-f -functions can be 
made to represent another -function with a required degree o-f 
accuracy. The term itsel-f o-f ten brings into mind the word 
“ perpendi cu 1 ar “ , and this thought can be a visualization o-f 
the structure o-f the members o-f an orthogonal set. 

Consider that we have a -function set, Spj<t), where n = 

9, 1,2, 3,... The set is said to be orthogonal with weight K 

over the interval 0>t>T i -f 

T K i -f n=m 

K-Sn<t)S^<t)dt = <B. 1) 

0 0 n^ 

with n and m being integer values. I-f the constant K is equal 

to 1 then the set is normalized and the set is re-ferred to as 

an orthonormal set. 

With this orthogonal set, we may now represent another 
function, f(t), defined over the interval <0,T), as 

T 

f<t) = SC^S^Ct) <B.2) 

0 

and is a number that indicates the value of the function 
Sf^(t). Cp, can be chosen so as to mimimize the mean-square 
error in representing f<t). 
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1/T 



(B.3) 



C 



n 



T 

•f < t) t) dt 
0 



This orthogonal function series representation reduces the 
number of coefficients needed to completely represent the 
signal. C23: pp.1-33 
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APPENDIX C 



MODULQ-2 ADDITION AND THE GRAY CODE 



Modulo-2 addition is an important mathematical operation 
in Walsh Theory. Its operation is used in the de-finition o-f 
the Pal ey ordered Walsh -functions with Rademacher -functions, 
and also in the product o-f two Walsh -functions, namely, 

WAL(n,t)WAL<m,t) =WAL(n«m,t> <C.l) 

where n'Sm indicates modulo-2 addition. 

Modulo-2 addition is really binary sums without the carry, 
and obey the rules 

0-H0=0, 0-U=l, U0=1, and 1-H=0 

The Gray Code is a binary code that is often used in 
communications because the codes for successive decimal digits 
differs by only I bit. It is not a weighted code, meaning 
that the decimal value of a coded digit cannot be computed by 
a simple formula. 

Table C.l. Gray Code for 16 Digits. 



Decimal 


Code 


Decimal 


Code 


0 


0000 


8 


1100 


1 


0001 


9 


110 1 


2 


0011 


10 


nil 


3 


0010 


1 1 


1110 


4 


0110 


12 


1010 


5 


0111 


13 


1011 


6 


0101 


14 


100 1 


7 


0100 


15 


1000 
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APPEhJDIX D 



FAST WALSH TRANSFORM AND POWER SPECTRAL DENSITY PLOTS 

The following pages are the plots of the FWT and 
PSD coefficient of the simulated TOA TAG strings used in the 
thesis. 

The plots are normalized to the maximum component and thus 
show the relative value of each component to this maximum. 

The value of the maximum component is indicated on the 
vertical axis. 
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Figure 0.1. FWT ; TCDA TAG 2. 
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Figure D.2. FWT : TOA TAG 3, 
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Figure D.3. FWT ; TOA TAG 4. 
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Figure 0.4. FWT : TOA TAG 5. 
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Figure D.5. FWT ; TOA TAG 6, 
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Figure D.6. FWT : TOA TAG 7. 
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Figure D.7. FMT : TOA TAG 8. 
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Figure D.3. FWT : TOA TAG 9 . 
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Figure D.9. FWT : TOA TAG 10. 
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Figure D.18. FWT ; TOA TA6 11. 
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Figure D.ll. FWT ; TOA TAG 12. 
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Figure 0.12. FWT ; TOA TAG 13. 
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Figure D. 13. FWT : TOA TAG 14. 
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Figure D.14. FWT : TOA TAG 15. 



189 



WALSH INDEX (n) 



FAST WALSH TRANSFM 
TOA TAG 16 




I 



_L. 

ID 



I 

1 

I 

♦ 






_1_ 

in 

f 

I 



(S290-T) DVW 13d 



(D 



(D 

in 



OD 






OJ 

CO 



c\J 

CO 



CD 



J 

— • 
I 



Figure D.15. FWT ; TOA TAG 16. 
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Figure D.16. FWT : TOA TAG 19. 
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Figure D.17. FWT : TOA TAG 22, 
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Figure 0.13. FWT : TOA TAG 25. 
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Figure D.19. FWT ; TOA TAG 32. 
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Figure D.20. FWT : TOA TAG 40. 
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Figure D.21. FWT ; TOA TAG 51. 
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Figure D.22. FWT : TOA TAG 64. 
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Figure D.24. PSD ; TCA TAG 3. 
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Figure D.26. PSD ; TOA TAG 5. 
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Figure D.27. PSD ; TOA TAG 6 . 
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Figure D.28. PSD ; TOA TAG 7. 
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Figure D.29. PSD : TOA TAG 8. 
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Figure D.38. PSD ; TOA TAG 9. 
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Figure D.31. PSD : TOA TAG 18. 
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Figure D.32. PSD ; TOA TAG 11. 
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Figure D.33. PSD : TOA TAG 12. 
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Figure D.34. PSD ; TOA TA6 13. 
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Figure D.35. PSD : TOA TAG 14. 
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Figure D.36. PSD : TOA TAG 15. 
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Figure 0.37. PSD : TOA TAG 16. 
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Figure D.38. PSD : TOA TAG 19. 
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Figure D.39. PSD : TOA TAG 22. 
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Figure D.48. PSD ; TOA TAG 25. 
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Figure D.41. PSD ; TOA TAG 32. 
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Figure D.42. PSD ; TC^ TAG 40. 
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Figure D.43. PSD ; TOA TAG 51. 



218 



SEQUENCY 



POWER SPEC DENSITY 
TOA TAG 64 



(88T00-T) 3VW 13y 



CM 

j 

J !2 
1 00 

-i CO 

I CM 

4 ^ 

i CM 

ts 

CM 

00 



CD 



- CM 

^ T— I 

4® 

i ” 

1 

■i 00 

L 

-} ^ 

-j 

' CM 



CJJ CD 



CD in 



cn CM 






Figure D.44. PSD ; TOA TAG 64. 
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Figure E.l. GROUP SPECTRUM : TOA TAG 5. 
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Figure E.3. GROUP SPECTRUM ; TOA TAG 14. 
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Figure E.4. GROUP SPECTRUM ; TOA TAG 19. 
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